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On the Adiabatic Couette Flow of a Compressible Fluid 


By Wiri1aM F. HuGues!) and J, FLETCHER OSTERLE?), Pittsburgh, Pa., USA 


Nomenclature 


P pressure (absolute) ; 

T temperature ° R; 

* distance along slider; 

L length of slider; 

fe viscosity ; 

C, specific heat at constant volume; 

M mass rate of flow; 

h distance between parallel plates; 

u fluid velocity in the x-direction; 

velocity of the slider pad in the x-direction; 
coordinate distance normal to the surface of the slider; 
density ; 

temperature coefficient of viscosity ; 
distance coefficient of viscosity; 

gas constant; 

specific heat at constant pressure; 

specific heat ratio. 


CP WRROEREG 


» inlet conditions at 7 = 0; 


k 
Subscripts: 
() 

(), outlet conditions at + = L. 


Introduction 


With the increasing use of air as a lubricant and interest in pneumatic 
control, a fundamental study of air flow in a narrow passage seems desirable. 
The problem investigated here is the adiabatic flow of air through the passage 
between two parallel plates very close together and in relative tangential 
motion, so-called Couette flow. The case where a pressure differential exists 
across the moving plate is also considered. 

In Couette flow, heat generated by the viscous shearing stresses in the fluid 
is partially stored in the fluid increasing its internal energy, and partly con- 
ducted away through the plates which bound the flow. The isothermal and 
adiabatic solutions describe then the two limiting flow conditions. Whether the 


flow is isothermal, adiabatic, or somewhere in between depends on the indivi- 


dual situation, but it is thought that in many practical situations the flow is 


1) Department of Mechanical Engineering, Carnegie Institute of Technology. 
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adiabatic [1]?). For isothermal Couette flow, a pressure gradient cannot occur in 
the film; it being possible only when the plates are tilted. However, for adia- 
batic flow, a pressure gradient is incurred for parallel operation of the plates 
because of the so-called ‘thermal wedge’ effect, which is a result of compres- 
sibility, and may become important for very high relative velocities [2]. _ 

The first treatment of the Couette flow of air was included in the work of 
HARRISON in 1913 [3] who considered the more general problem of isothermal 
air flow between tilted plates. Work was extended to a consideration of momen- 
tum effects in air flow for stationary plates along with experimental data by 
S. GRINNEL in 1955 [4]. The following mathematical analysis allows the pre- 
diction of the pressure and temperature distributions and mass rate of flow for 
adiabatic Couette flow, and certain conclusions are drawn as to the importance 
of the thermal wedge effect. 


Mathematical Analysis 


For the viscous laminar flow of a fluid through a narrow channel! the usual 
assumptions common to lubrication type flow can be made [1]. In addition, for 
air, the compressibility terms are neglected in the momentum equation; this 
assumption is valid as long as the velocities are much less than sonic [5]. The 
momentum equation then takes the simplified form: 

2 
a = bh sot (1) 


The energy equation for air may be written [6]: 
(2) 


Since the viscosity of air varies nearly linearly with temperature through a 
range of several hundred degrees, the following relationship is assumed [7]: 


f= oll + y (LT — T)] . (3) 
The equation of state f 
PHoRT (4) 


ge T a Pau (Mh 


must be used, and the continuity equation written in the form 


h 
M = [oudz=const (5) 
0 


will be employed. 
Equations (1) through (5) must be solved simultaneously to completely 
describe the flow. 


*) Numbers in brackets refer to References, page 95. 
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By integrating equation (1) subject to the boundary conditions: 


which arise from the ‘no-slip’ assumption of hydrodynamics, we obtain 


A etic 25 18 Uz 
“4 = X Fe (22—hz)+ oe (6) 


Zu 


for the velocity profile. Substituting this expression and the equation of state 
into the continuity equation there results 


he dP Uh MRT 
12 pu A aa 3 B ae _ (7) 


Now, substituting equation (6) into the energy equation and integrating we 
obtain: 


ve a2; Oo (<\" wU? 


© dx 12p \ds. fot 


mS d i or (8) 


Raa Dah we ipiuige 
The three equations (5), (7), and (8) constitute a set of equations to be solved 


for the temperature, pressure, and viscosity variations. 
Eliminating T between equations (7) and (8), making use of the relationships: 


Ge tek | oreo ie 
Cy 
there results 


k h® wae (= oc) 16 Us" Uh dP 


12(k—1) * dx \p * dx h 2k =i) “ae (9) 


If the viscosity variation with x were known equation (9) could be integrated 
directly for the pressure distribution. If a linear variation of viscosity with x 
is assumed, it will be found that the final temperature distribution is nearly 


Bearing 


Figure 1 
The parallel-surface slider pad and bearing in Couette flow. 
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linear in « and, hence, the assumption is a reasonable one for a first approxi- 
mation. We assume 
[b= fg (1+ a x) , (10) 


where « is determined by the maximum temperature at x = L. 
Equation (9) can now be integrated directly to yield 


adP My (1 + & x) My (1 + & #) ax\ | 
Sige: al Ra SF Ae et ee B-o% (1 = “74) [? (11) 
where 
a WN IU (ee ah), 
a 
and 
k h? 
Com stig 


and B is a constant of integration to be determined by the boundary conditions. 
In order to integrate (11), let 


= w(x) 6 f(x), (12) 


where 


f(x) = Bo (x+ 22). (13) 
Then we obtain for w(x) ) 
dw w 6? De be aes 
dx ieee T(*) , ty) 


The integral of equation (14) has three forms depending on whether k (k — 1) 
is greater than, less than, or equal to unity. For air, k (k — 1) is eee to 
0-56, and if we let 


4 


2 
a it 5 


kh (& —1) 


equation (14) integrates to the form 


(1 — a?) w+ 2(1— a) 
(1 — a?) w+ 2(1+ a) 


in= yi — 2) w?+4w+ 4] (f(a)? = 2 In | 


|; 05) 
where C is a constant of integration. The boundary conditions are the following . 
tO; Ps Pye le ate alo ae 


Applying these conditions to equation (12), the boundary conditions on w(x) are 


iP 12 
pee = : aN fae ain Cone 
x w BO: Doe wm ors) Of(L)’ 


These conditions may be applied to equation (15) to obtain two equations 
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involving the constants B, C, and «. 


(t—a@*) (P,\2 , PoB vate 1 (= @*\ Pi + 2 at 
in gon fen + C29 + ad ae om |: — a’ 2B (i Bd (16) 


{[G= at) Pe, Pe HZ) HE) _ 3 (1 — a*) P, + 26 f(L) (1 —a) 
I jeaaenes Aas Pare’ be —a¥) P, + 20f(L) (+ ay | 


We can obtain another equation relating B and « by combining equation (7) 
with equation (11) to eliminate dP/dx. Hence 
re) 


eee ' Uh— us [P + 0 f(x] (18) 


From the x and T viscosity variations it follows that 
T=——+Y,. (19) 


Now, by combining equations (18) and (19) and using the boundary conditions 
on pressure and temperature it follows that the mass rate of flow is given by 


Uh h? (1 B 
N= 00 [- 2 12 c 3 al ae 
and 
24 P, es ne =| Uh oC eh 
Bae Ley, - 2 I2\.0 rs Bee, UENO P ! (at) 
{1 24 PR Qo oe wie Eat 
vor®L | 2 HW h ie le = | 


For the special case of Py, = P,, that is the non-pressurized slider pad,.equation 
(21) takes the simpler form: 


BAR By f Udy ht B \=2 a 
Saabe gael 2 alison ag | =f : a 


We can eliminate C by combining (16) and (17) to obtain 
[(1 — a?) P? + 4B6P, + 4.B26?) B? 
{(1— 4%) P2 + 4 P,0 f(L) + 46? (HEN*} 
1 (fl 2) Py + 2B 0 alli —a*) P) 220 ALY + a) 


=| "aa P, + 2B Oi +a) [i — a4) P, +26 7(L) (1 — a) * 


In 


The use of these equations for the solution of practical flow problems is straight- 
forward, but involves some trial and error work since the equations are nearly 
all transcendental. Equations (24) and (21), or (22), can be solved simultane- 
ously for B and «; however, since the unknowns do not appear in explicit form, | 
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a trial and error numerical solution must be effected. The constant C may then 
be determined explicitly from equation (16) or (17). Once «, B and C are known, 
equation (15) gives the pressure distribution and equation (18) gives the tem- 
perature distribution in terms of the pressure and mass rate of flow which 
may be determined from equation (20). 

It is interesting to note that for the hydrostatic situation of zero relative 
velocity of the slider pad there is no increase in temperature of the air as it 
flows through the passage. This may be shown by examining equation (9) in 
the limiting case of U = 0: 


k h8 Gaal dP 
ates ered eiae eipats “a 
Therefore 
ei x cee const . (25) 
bh ax 
‘And from equation (7) 
J AS 
ax ; 


From these equations it is seen that the flow must be isothermal, consequently 
isoviscous, and by integrating equation (7) under these conditions the pressure 


distribution is found to be 
la ee cen 2 
Pepe sf 


(26) 
which is that obtained by GRINNELL [4] for the isothermal case. 


Numerical Example and Discussion 


The analysis herein presented allows the prediction of pressure and tem- 
perature distributions and the mass rate of flow for adiabatic Couette flow. 
To illustrate the use of these equations sample calculations were made for the 
' following slider pad: 


Ee = nee 
fo 0-001 ins 
U = 1000 in./s; 
dg = 560" KR; 
i, = 147 psia, 
P= 14-7 psia, 


y =1:33x10-3/°R 


and the results for temperature and pressure distributions are shown in Figure 2. 
The temperature which undergoes appreciable rise, 146° R, is nearly linear in 
x, and the pressure undergoes a slight build up with a maximum rise of about 
1:9 psi due to the thermal wedge effect. Hence, if a pressure differential is 
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maintained across the moving plate, or slider pad, the contribution to the pres- 
sure by the thermal wedge effect may be quite small compared to hydrostatic 
pressures. This is the case in hydrostatic lubrication where hydrodynamical 
effects may be negligible for low speed operation. 

For very high speed operation the thermal wedge effect becomes important 
and may be used to provide a load supporting pressure build up in a journal 


Pressure (P-F,) psi. —= 
Temperature (T-Tg) °F. —=— 


x Distance [in.}) —=— 


Figure 2 
Pressure and temperature distributions. 


or slider pad bearing in addition to geometrical wedge effects which are present 
for isothermal operation of tilted pad slider bearings. 

In actuality the operating characteristics for a slider pad bearing falls 
somewhere between the isothermal and adiabatic limits, depending on the 
conductivity of the bearing material and degree of insulation. For most prac- 
tical situations, the flow will be nearly isothermal because of the relatively 
high conductivity of the bearing material. Recent work by the authors seem 
to indicate in a quantitative manner that the flow may be essentially isothermal 
in many cases, but the adiabatic solution does form an upper bound which 
might be of considerable interest in certain problems. 
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Zusammenfassung 


Eine Untersuchung wurde durchgefiihrt tiber die adiabatische Stromung der 
Luft zwischen zwei parallelen Platten, bei einem Druckunterschied quer tiber die 
sich bewegende Platte. Die Losung lasst eine thermale Keilwirkung der hydro- 
dynamischen Schmierung erwarten und gibt eine quantitative Analyse des 
allgemeineren Problems, wo ein Druckunterschied quer iiber die Gleitplatte vor- 
handen ist. Die Ergebnisse deuten an, dass der Beitrag des hydrodynamischen 
Druckes oft dusserst klein sein kann (fiir massige Geschwindigkeiten) im Ver- 
gleich zum hydrodynamischen Druck in den dem Druck ausgesetzten Lagern. 
Im Falle von stillstehenden Platten ist die adiabatische Lésung gleich der iso- 
thermen Loésung; weswegen keine TemperaturerhGhung stattfindet. 


(Received: September 24, 1956.) 


Comment agit le cable tracteur d’un téléphétique 


Par Marie Marcet Nicoras, Paris?) 


1. Description générale schématique d’un téléphérique 


Un téléphérique se compose essentiellement d’un cable tendu entre deux 
points A et B sur lequel roule un chariot auquel est suspendue une cabine. Le 
déplacement du chariot sur le cable porteur, véritable rail aérien, est obtenu 
par l’action d’un cable tracteur. 

Le présent mémoire a pour but l’étude de l’action d’un cable tracteur pro- 
voquant le déplacement d’un chariot de téléphérique sur un cable porteur. De 
ce dernier nous ne rappellerons qu’une caractéristique mécanique fondamen- 
tale; a l’une de ses extrémités le cable porteur est soumis a une traction cons- 
tante. De ce fait toutes les caténes que peut former ce cable entre les appuis 


1) Secrétariat d’Etat aux Travaux Publics, aux Transports et au Tourisme. 
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A et B, et éventuellement entre ceux intermédiaires, sont toutes découpées sur 
des chainettes de modules pouvant étre différents mais dont toutes les bases 
canoniques sont a la méme altitude, ce que nous exprimons en écrivant que 
toutes les caténes de suspension sont isosthénes (en un point quelconque de 
lune d’entre elles la traction est mesurée par le poids de cable ayant une 
longueur égale a l’altitude du point au-dessus de la base canonique commune). 

Comme on le verra par la suite, cette particularité mécanique du cable 
porteur a une importance considérable dans l'étude du cable tracteur qui lui 
est associé. 

Les téléphériques 4 va-et-vient sont ceux dont le chariot parcourt un cable 
porteur tantdét dans un sens de A vers B, tantét dans le sens opposé. Il est donc 
normal de mettre en paralléle deux cables porteurs et de faire en sorte qu’avec 
un seul cable tracteur on assure simultanément les déplacements des deux 
cabines en sens opposes. 

Il existe des téléphériques a va-et-vient dans lesquels la deuxiéme voie 
paralléle n’a pas été installée; méme dans ces installations l’aménagement et le 
fonctionnement du cable tracteur demeurent analogues a ceux a deux chariots, 
que nous décrivons ci-aprés. 

Le cable tracteur se compose d’un brin fixé a la face cété B du chariot; ce 
brin va contourner une molette (ou une poulie) dans la station B pour revenir 
a la face coté B de l'autre chariot a laquelle il est fixé; de méme un autre brin 
issu de la face c6été A du premier chariot va contourner une poulie (ou une 
molette) dans la station A pour atteindre la méme face du second chariot ot 
il est fixe. 

Ces deux brins constituent donc un cable continu formant une boucle entre 
les extrémités A et B, avec les deux maillons intercalés que sont les chariots. 

Il n’est pas exclu que le brin issu par exemple de la face B du chariot ne 

‘rejoigne pas celle de l’autre chariot immédiatement aprés avoir contourné la 
molette (ou la poulie) B; il existe en effet des installations ot un tel brin ne 
rejoint la face B de l’autre chariot qu’aprés avoir contourné une poulie (ou 
molette) en A, puis une deuxiéme molette (ou poulie) en B (différente de la 
premiére); ce dispositif, di 4 M. l’ingénieur REBUFFEL, permet d’utiliser les 
parties BA et AB de ce circuit comme cable d’agrippement respectif du second 
et du premier chariot. 

On concoit qu’avec un tel aménagement une molette motrice installée dans 
l'une quelconque des deux stations A ou B, et autour de laquelle le brin du 
cable tracteur sera convenablement enroulé, suffira pour assurer les déplace- 
ments du chariot. Pour limiter les eftorts de traction auxquels sera soumis le 
‘cable tracteur on organise dans l’une des stations un dispositif de mise en 
tension au moyen d’un contrepoids. 

| Ce dispositif de mise en tension peut étre organisé indifferemment dans la 

station ou se trouve la molette motrice, ou dans celle opposée dont la molette 


_ZAMP VIII/7 
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n’est qu'une poulie de renvoi, folle autour de son axe. Dans ce dernier cas la 
mise en tension s’effectue en exergant une traction constante, par le moyen 
d’un contrepoids, sur l’axe qui peut se déplacer entre deux glissiéres. L’effort 
de traction constant exercé sur le brin du cable dans ce systéme est donc égal 
Me ad en 

Dans l'autre cas, en appelant A la station ot sont réunis la mise en mouve- 
ment et la mise en tension, le brin tracteur, issu de la face A du chariot, en 
pénétrant dans cette station, A, contourne d’abord une poulie montée folle 
autour d’un axe horizontal (supportant un contrepoids P), puis contourne la 
molette motrice, oti il recoit une impulsion, et revient contourner une deuxiéme 
poulie montée folle, jumelle de la précédente mais non solidaire de celle-ci, 
autour de l’axe support du contrepoids; aprés ce circuit, tracteur-tenseur, dans 
la station mixte A, le brin tracteur va rejoindre la face A de l’autre chariot. 
Dans ce systéme les efforts de traction subis par le brin tracteur aux points ot 
il entre dans la station A et ou il en sort, ne sont pas égaux, l’un étant égal a 


Vautre a 


en appelant M l’effort transmis au brin tracteur par la molette motrice, lequel 
effort peut s’exercer indifféremment dans le sens du défilement du brin tracteur 
ou dans le sens opposé. 

Si donc on connait l’effort maximum M que la molette motrice pourra trans- 
mettre au brin qui la contourne, on en déduit la valeur qu’il faut donner au 
poids du contrepoids pour que la traction dans le cable tracteur dans la station 
mixte demeure inférieure a une limite donnée. 

On remarquera que cette installation mixte comporte une imperfection 
mécanique. En effet, les deux poulies folles jumelles exercent sur l’axe hori- 
zontal support du contrepoids P deux actions verticales qui ne sont pas égales 
lune étant 


l’autre étant 


Il en résulte une action de basculement de cet axe, mesurée par le moment 
M e¢ ote est l’écart entre les deux plans verticaux axiaux des poulies jumelles. 
Au cours d’un cycle de parcours, il est certain que les tractions exercées dans 
les brins tracteurs s’échangeront, donc la sollicitation au basculement changera 
de sens, ce qui nécessiterait un aménagement spécial, pratiquement irréalisable, 
des glissiéres guides de cet axe d’une horizontalité approximative variable. On 
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peut douter que cette imperfection mécanique des stations mixtes, tenseur- 
moteur, soit favorable 4 un défilement stable de leur cAble tracteur. 

En conclusion de cette description de l’aménagement général du cable 
tracteur d’un téléphérique on retiendra que celui-ci se compose de deux brins 
distincts mis bout a bout, avec interposition de deux maillons constitués par 
les chariots, pour constituer une boucle continue entre les stations A et B. 


Molette motrice supportée par un axe 
reposant sur un bati fixe 


Poulies jumelles de renvoi dont l’axe 
commun soutient le contrepoids 


Contrepoids 


Figure 1 


Schéma du dispositif moteur—tenseur. 


Celui de ces deux brins dont les deux extrémités sont fixés a la face tournée 
vers la station B, de chacun des chariots sera dit: brin B; et de méme l'autre 
brin sera dit: brin A. Un seul de ces brins contourne la molette motrice, l’autre 
ne contournant qu’une poulie folle: le premier sera qualifié brin actif, et le 
second brin rétroactif. 

Celui de ces deux brins, qui contournera la ou les poulies de mise en tension, 
cumulera donc son role de brin, actif ou rétroactif, avec celui de brin de mise 
en tension. 

Il y a donc deux combinaisons possibles pour associer la mise en tension 
soit avec la mise en traction soit avec la mise en rétrotraction. Et de fait on 
trouve dans la pratique ces deux combinaisons, sans égard au sens de la diffé- 
rence d’altitude entre les deux stations. 

On remarquera que pour décrire ]’aménagement général d’un téléphérique, 
et plus particuliérement celui des cables tracteurs, nous avons évité d’évoquer 
les différences d’altitude entre les deux stations, ce qui aurait été le cas en 
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qualifiant l’une de supérieure, l’autre d’inférieure. En effet cette différence 
d’altitude n’est nullement une caractéristique sine qua non d’une installation 
téléphérique, et nous soutenons qu’est encore une installation téléphérique 
celle destinée par exemple au franchissement d’un large ravin moyennant un 
cable porteur tendu entre deux points au méme niveau. 

De méme dans la succession des appuis intermédiaires, il peut fort bien 
s’en trouver trois successifs, dont celui du milieu soit le plus élevé. 

Dans les deux cas susvisés une distinction procédant du critérium d’une 
moindre sollicitation d’un des deux brins ne peut s’appliquer. 


2. Définition du probleme mécanique du cable tracteur 


Ces préliminaires étant posés, nous pouvons aborder l|’étude du compor- 
tement du cable tracteur, plus exactement de ses deux brins quant a leur 
intervention dans le déplacement en va-et-vient, des deux chariots. Si ces 
interventions n’étaient pas comparables, on pourrait déja conclure qu’en raison 
de la symétrie inversée des phases successives, au cours d’un cycle de parcours, 
chaque brin assurerait tour 4 tour un role tracteur prépondérant. 

Il n’en est rien. Comme nous l’avons écrit dans une étude antérieure (Cable 
de téléphérique supportant une charge, Ann. Ponts Chaussées 7955, mai—juin) 
il est injustifié de minimiser en apparence l’utilité d’un brin du cable tracteur 
en lui attribuant la désignation de cable-lest a l’opposé de celle de cable 
tracteur réservée a l'autre. 

Depuis que nous avons dénoncée Vimproprieté de cette dénomination 
«cable-lest» nous en avons vainement cherché une définition précise. Par 
exemple dans le cas d’un téléphérique célébre, récemment mis en service, ou 
la station a l’altitude inférieure est une station mixte moteur-tenseur, il est 
évidemment impossible de soutenir que l’un des brins du cable tracteur ait un 
role prépondérant de tracteur, tandis que l’autre, qui deviendrait dénommé a 
tort cable-lest, ne serait qu’un simple ballot servant seulement a régler la 
traction du brin affublé du titre de cable tracteur. 

L’étude de l’effort qu'il est nécessaire d’exercer sur un chariot roulant sur 
un cable porteur pour provoquer son déplacement montre: 

le lVimpossibilité d’avoir la maitrise de ce déplacement au moyen d’un 
cable tracteur n’agissant que par un seul brin; 

2° la rigoureuse équivalence des réles tracteurs de chacun des deux brins, 
en tout point du parcours du chariot. 

Par «avoir la maitrise du déplacement d’un chariot roulant sur un cable 
porteur» nous entendons avoir la possibilité de choisir, a l’intérieur de certaines 
limites, la trajectoire de ce chariot et celle de lui imprimer sur cette trajectoire 
telle vitesse et telle accélération, toujours dans certaines limites, dont en parti- 
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culier une vitesse nulle, c’est a dire maintenir le -véhicule suspendu a l’arrét en 
un point quelconque de son parcours. 

Les limites du choix de la trajectoire sont celles dont nous avons démontré 
lexistence (Ann. Ponts Chaussées 7955, mai-juin, 338). Tout se passe en effet 
comme si le chariot, réduit 4 un point matériel, était astreint 4 se mouvoir dans 
un champ de forces résultant de deux autres champs de forces, savoir: 1° les 
forces verticales de pesanteur agissant sur une masse égale a celle du chariot et 


(BS eo 


//// Zone A une force applicable par point. 
//## Zone a deux forces applicables par point. 


Figure 2 


de sa charge; 2° les forces que le cable porteur déformé en deux caténes AC et 


BC exerce sur ses deux sections qui encadrent le point C. Ce deuxiéme champ 
de forces indépendantes de la masse du point est déterminé par la longueur 


et l’inclinaison de la sustente AB, par le poids par unité de longueur du cable 
porteur et enfin par la traction constante (contrepoids) exercée a lune des 
extrémités. 

Nous reproduisons le dessin qui résume cette partie de cette étude; on 
remarquera qu’y est indiquée la direction générale de la force que le chariot 
devra exercer sur le cable porteur pour contraindre celui-ci a se déformer en 
deux caténes dont les extrémités communes marquent l’emplacement du 
chariot. 

Cette force bien définie en chaque point de l’espace accessible au chariot ne 
peut provenir que du poids de celui-ci et de la transmission, par le chariot, de 
Yeffort qu’il subit de la part du cable tracteur. Donc en définitive l’effort que 
devra exercer le cable tracteur sur le chariot pour maintenir celui-ci a l’arrét 
en un point quelconque de la zone d’accessibilité est rigoureusement déterminé. 
Pour un méme point d’arrét, mais pour des poids variables du chariot Veffort 
du cable tracteur gardera une composante horizontale invariable (égale a la 
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différence des modules des deux caténes tracées par le cable porteur déforme) ; 
seule variera la composante verticale de cet effort. 

Si de plus on se donne la trajectoire et la vitesse (variable ou non) le long de 
celle-ci, effort demandé au cable tracteur sera encore entiérement déterminé 
compte tenu de la vitesse, du rayon de courbure de la trajectoire et de l’accélé- 
ration tangentielle. 


Figure 3 
Couple de caténes de suspension appareillées pour une sustente 4B. 
bs ba livet R a a HN Boa, Fee 
= ‘ = A — eae et eee 
A, =~ X Aa; By =p X Bb, Nabe gaa On Oh ie area 


En définitive le probléme mécanique que nous avons a résoudre est le sui- 
vant: étant donné un point C de l’espace, comment peut-on avec un cable 
tracteur issu de deux points fixes A et B le faire passer par un troisiéme point C, 
pour qu’il y exerce une force dont les deux composantes, horizontale H et verti- 
cale V, aient des valeurs imposées. 

La solution de ce probléme nous sera donnée par |’étude des tractions 
diverses qu’un cable, astreint 4 passer en deux points fixes, exerce sur ces 
points selon les diverses caténes de suspension qu’il peut former. 


3. Etude des tractions exercées par un cable fixé en deux points 
A et B entre lesquels il est suspendu 


On sait que toutes les caténes de suspension possibles entre deux points A 
et B peuvent étre associées deux a deux, l’une tendue de module 7;,, l’autre 
lache de module 7; (7; <7;), et ayant toutes deux méme base canonique. 


Cela signifie en particulier que les ordonnées canoniques A a, et Bb, des 
extrémités ont mémes valeurs dans les deux systémes de coordonnées canoniques 
afférentes aux deux chainettes. Les abscisses canoniques de ces extrémités 
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étant respectivement O” ay et O” Dy en suspension lache et O’ ie O' by en 
suspension tendue. 


On sait aussi qu ‘il existe une suspension dans laquelle les ordonnées cano- 


niques A a et Bb sont minimales et que cette suspension de module 7 sépare 
la zone des caténes tendues de celles des caténes laches. 


Comme A a, Ax, B b,., etc., mesurent a un coefficient constant pres (le poids p 


du cable par unité de longueur) les tractions Ax et B,, exercées par le cable sur 
les extrémités, on justifie ainsi la dénomination de suspension minimale pour 


celle de module y. Rappelons que la suspension minimale d’une sustente AB 
faisant un angle # avec l’horizontale est définie par la condition que les tangentes 
a la chainette en A et B se coupent sur la base. Analytiquement cette condition 
se traduit pas les deux relations suivantes entre les abscisses canoniques Xp, %; 
des extrémités, dans la chainette de module 7: 


(voir a ce sujet notre mémoire Equilibre d’un cable de téléphérique, Ann. Ponts 
Chaussées 7954, janvier—février). 


Les tractions A, comme les tractions By, sont égales mais avec des inclinai- 
sons différentes, celles afférentes a la suspension lache étant plus voisines de la 
direction verticale. 

Convenons d’affecter la lettre B a Vextremite supérieure, de désigner par 


V;, la composante verticale de la traction A, en suspension tendue et par 
Vj, la méme composante en suspension lache; de méme nous désignerons par 


W;. et Wz les composantes verticales de la traction B, en suspension respecti- 
vement tendue et lache. Les composantes horizontales H;, des tractions A, et 


B,, en suspension tendue étant égales a 7;,, mais dirigées en sens contraires 
(toujours vers l’aplomb de la sustente); il en est de méme des composantes 


horizontales Hy (égales a r;) des tractions A, et B, en suspension lache. 

Une étude détaillée, que l’on trouvera en appendice au présent mémoire, 
apprend que des relations, d’une part entre V et H, d’autre part entre W et H 
il existe des représentations graphiques dont les caractéristiques essentielles 
suffisent pour répondre le probléme mécanique du cable tracteur tel que nous 
Vavons posé. 

En effet si d’un point quelconque « (pour lequel et opportun de choisir 
l’extrémité inférieure A) on trace tous les vecteurs A, figurant les tractions 

_exercées par les diverses suspensions sur ]’extrémité inférieure A, leurs extre- 
mités ponctuent une courbe dont la concavité continue est orientée vers la 


104 Marie Marcer NICOLAS ZAMP 


zone inférieure sous l’aplomb de B. Les abscisses de cette courbe (A) mesurent 
en chaque point le module de la suspension (ou la composante horizontale Hy, 
avec son orientation), tandis que les ordonnées figurent en grandeur et sens 
les composantes verticales. Nous appelons le point « origine des vecteurs A, le 
foyer de la courbe (A). 


Figure 4a 


(i = 28° 29’ 30”). Représentation des courbes de traction (A) et (B) aux extrémités d’une sustente 
d’inclinaison 7 et de projection horizontale A by. Deux verticales équidistantes du milieu de linter- 


valle A by déterminent les tractions Aet Be d’une méme caténe de suspension. Les vecteurs de traction 
sont figurés avec une longueur égale a la moitié de la longueur de cable ayant méme poids. 

Nota: Sion tragait les deux courbes a partir d’une méme verticale et en retournant par exemple 
lune d’elle (B), la courbe (—B) n’aurait aucun point d’intersection avec la courbe (A) étant. partout 
au-dessous de celle-ci. 


Cette courbe concave a deux asymptotes, l’une verticale (orientée vers le 
bas) pour 7 = 0, c’est a dire pour une caténe de suspension infiniment lache, 
l'autre faisant au-dessus de Vhorizontale un angle 7 égal a celui de la sustente, 
c’est a dire pour une caténe de suspension infiniment tendue, le point d’inter- 
section de ces deux asymptotes ayant une abscisse nulle et une ordonnée égale 
au vecteur-force, figurant un poids de cable d’une longueur égale a la moitié 


de celle I’ de la sustante AB. 
Cette courbe comporte deux points remarquables, savoir celui A,, ow elle est 


normale au vecteur re du point considéré, et celui A, ot. la composante verticale 
s'annule et change d’orientation. Ces deux points peuvent, pour chaque incli- 
naison 7, étre déterminés sans difficulté et avec précision puisqu’ils correspon- 
dent, le premier a la traction en suspension minimale, le second a la traction 
en suspension dont la tangente a l’extrémité inférieure est horizontale. 
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Si a partir d’un autre point , pour lequel il est opportun de choisir le point 
by 2 aplomb de l’extrémité B sur Vhorizontale de l’extrémité inférieure A, nous 
déterminons de méme la courbe tracée par les extrémités B, des diverses trac- 
tions, nous obtenons une courbe (B) 4 concavité continue orientée vers la zone 


inférieure sous l’aplomb de A et comprise entre deux asymptotes ayant mémes 
dispositions que celles indiquées pour la courbe (A). Nous appellerons le point 


origine des vecteurs B,, le foyer de la courbe (B). 

La traction Bn de la suspension minimale marque le point B,, ou cette 
courbe est normale au vecteur Be. générateur de ce point. 

La détermination du point A, (avec une composante horizontale H, = A.) 
aura permis de situer sur la courbe (B) le point dont le vecteur B, a méme 
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composante horizontale (méme module) H,, mais orientée en sens inverse. On 
démontre qu’entre ce point B,, et ce point B, la courbe (B) a-une tangente hori- 
zontale: ceci signifie que c’est pour une suspension comprise entre la suspension 
minimale et la suspension a tangente inférieure horizontale que la composante 
verticale W a une valeur minimale (donc en suspension tendue et a tangente 
inférieure plongeante). 

Il est évident que si la sustente AB change de longueur, mais en gardant 
la méme inclinaison 7, les deux courbes susdécrites (A) et (B) se modifient 
homothétiquement, dans le rapport des longueurs des sustentes a partir des 
origines de leurs vecteurs générateurs. 

Plus intéressante pour le probléme étudié est la déformation progressive des 
deux courbes (A) et (B) quand l’angle # varie, la projection horizontale de la 
sustente AB demeurant inchangée. 

On démontre alors que quand 7 augmente la nouvelle courbe (A) est au- 
dessus de la précédente, et qu’en particulier le vecteur he qui marque la sus- 
pension minimale (ou le point ot la courbe (A) est normale au vecteur A,,) 
diminue de longueur en méme temps que diminue l’angle qu'il fait avec l’hori- 
zontale, ces deux quantités tendant simultanément vers zéro, la forme de courbe 
(A) tendant vers celle d’une droite verticale. La longueur maximale de ce 
vecteur minimal A, étant celle afférente 4 une inclinaison 7 nulle, l’angle du 
vecteur A, avec l’horizontale étant alors égal a 56°28’ et A, étant égale au 
facteur p, prés a 0,726 fois la longueur de cette sustente horizontale. Mais 
pour 7 nul, la courbe (B) est devenue identique, mais inversée, par rapport a 
la verticale, a la courbe (A). En effet pour des valeurs croissantes de ?, les 
courbes (B) successives se placent les unes en dessous des autres entre deux 
branches asymptotiques dont l'une est toujours verticale tandis que l’autre 
s’incline en tendant a se rapprocher de la verticale. Et on démontre que le 
point haut de chacune de ces courbes ne cesse de s’abaisser en tendant vers 
le point inférieur a V’infini sur la verticale. (Pour 7 = 0 ce point haut se trouve 
sur l’horizontale d’ordonnée A 0/2 et a l’infini sur cette horizontale.) 

Les deux figures 4a et 4b représentent pour deux inclinaisons différentes, 
mais pour une méme longueur de la projection horizontale A by de la sustente, 
les courbes (A) et (B). 

Pour étudier l’effort de traction qu’un brin tracteur issu de A peut exercer 
sur un point C, nous tracerons la courbe (C4) afférente a l’extrémité C de la 
sustente fictive AC; de méme pour étudier l’effort de traction qu’un brin 
tracteur issu de B peut exercer sur ce méme point C, nous tracerons la courbe 
(Cz) afférente a l’extrémité C de la sustente fictive CB. 

Quel que soit le brin tracteur considéré on voit immédiatement que s'il 
agissait seul il ne pourrait exercer en C que trés exceptionnellement une force 
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imposée a priori. En effet seraient déja a éliminer toutes les forces imposées 
dont la composante horizontale serait orientée vers l’extérieur de l’'aplomb de 
la sustente fictive du brin tracteur. A la composante horizontale, convenable- 
ment orientée, et ayant une valeur imposée, il ne correspondrait dans le brin 
tracteur qu'une seule composante vérticale V. 

Il est évident que de telles conditions ne sont pas celles d’une maitrise, si 
limitée soit-elle, du mouvement d’un téléphérique. 

Supposons que par hasard cette composante horizontale de la traction du 
brin tracteur soit égale et opposée a celle de la force créée en C par le cable 
porteur déformé, supposons encore que la composante verticale de cette der- 
miére force associée au poids du chariot soit égale et opposée a la composante 
verticale V de la traction, on pourrait alors tenir l’arrét en ce point C avec le 
brin tracteur. Mais il suffirait de modifier le poids (la charge suspendue) du 
chariot pour étre assuré que le point C ne pourra plus étre un point d’arrét. 

Au contraire si on dispose de deux brins tracteurs, issus de points A et B, 
dont les verticales encadrent un point C accessible au chariot, on peut toujours 
réaliser théoriquement en ce point n’importe quelle force, le choix de celle-ci 
n’étant restreint que par une limite a la résistance des cables. C’est ce que nous 
allons démontrer. 


4. Réalisation d’une force quelconque au moyen de deux brins tracteurs 


Nous continuons a4 désigner par les mémes lettres, A et B, les appuis du 
cable porteur et les extrémités des brins tracteurs; en fait ces deux catégories 
d’organes mécaniques de I’installation ne coincident pas et ne peuvent coincider. 
Mais cette distinction n’est nécessaire que pour étudier le probléme plus parti- 
culier de l’abord ou du départ d’un appui par le chariot, aussi ne la ferons-nous 
pas dans l’exposé qui va suivre. 


Tout point C accessible sous une sustente AB (A plus bas que B) faisant 
avec l’horizontale un angle 7, marque le sommet d’un contour polygonal ACB, 


dont le cété AC fait avec l’horizontale un angle 7, toujours inférieur a /, si le 
point C est a une altitude intermédiaire entre celle de A et de B, tandis que le 


coté CB fait avec Vhorizontale un angle 4, toujours supérieur a 7 donc a 4,. 
En d’autres termes cette disposition des sustentes fictives signifie que C est 
extrémité supérieure pour la sustente fictive issue de l’appui inférieur et est 
extrémité inférieure pour la sustente fictive issue de l’appui supérieur. 

Mais il peut advenir, et ce sera presque toujours le cas si la différence 
d’altitude entre B et A n’est pas trés accentuée, que l’on ait a étudier des par- 
cours dont des points soient a une altitude inférieure a celle de l’appui le plus 


bas; dans ce cas la sustente fictive AC a une inclinaison 7, négative. Une telle 
disposition des sustentes fictives, caractéris¢e par un angle 7, négatif et un angle 
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i, positif, ou par le fait que C est extrémité inférieure pour les deux sustentes 
fictives, se rencontrera toujours si A et B sont a la méme altitude. 

Nous allons examiner dans ces deux cas comment le concours des deux 
brins tracteurs permet de réaliser en C n’importe quelle force donnée a priori. 


1° Cas d’un point C a une altitude intermédiaire entre les deux appuis A et B 


Nous mesurerons positivement les composantes horizontales des forces quand 
celles-ci sont orientées dans le sens aplomb de A (appui le moins haut) vers 
aplomb de B (appui le plus haut). 

Nous appellerons H, et W, les composantes, horizontale et verticale, de la 
traction CG exercée en C par le brin tracteur issu de A; nous appellerons Hz 


et Vz les mémes composantes de la traction Ce due au brin issu de B. Les 
composantes W et V sont mesurées positivement quand elles sont orientées 
vers le bas. 

A partir d’un foyer y, que pour la commodité de l’interprétation des courbes 
il y a intérét a choisir sur la verticale du point C étudié, nous tragons la courbe 
(C4) dont les coordonnées seront H4 et W, et la courbe (Cz) dont les coordon- 
nées seront Hz et Vz. 

Ces deux courbes sont en quelque sorte dos a dos de part et d’autre de 
Vaxe vertical des ordonnées V (ou W), lequel axe est aussi leur asymptote 
commune (dans ce cas les caténes tractrices ont un module nul, c’est a dire 
sont des caténes infiniment laches). Les autres asymptotes sont pour la courbe 
(C4) la droite faisant l’angle 7, avec V’horizontale et passant par le point de 
coordonnées: 0 et /'/2; pour la courbe (Cz) ce sera la droite faisant langle 7, 
et passant par le point de coordonnées: 0 et J’,/2. 

Dans la suite de notre exposé nous devrons considérer des courbes (C4) et 
(Cz) dont les asymptotes verticales ne sont plus confondues, aussi dans ce cas 
nous désignerons par y4 le foyer de la premiére et par yz celui de la seconde. 

Soient alors H et V les deux composantes d de la force ¢ a réaliser en conjuguant 


les deux tractions C, ace a: 3on devra avoir: H, H4 + H, B= =H. 


De méme on devra avoir: Wi + Vee an 

On réalisera une représentation graphique de ces deux relations: 

a) en déplacant de H (en grandeur et signe) la courbe (C4) puis en inversant 
celle-ci autour de sa nouvelle verticale asymptote, ce qui nous donnera la courbe © 
(—C,) inversée (le foyer y, sera déplacé de H par rapport au foyer y,) ; 

b) en tragant ensuite la courbe symétrique de (Cg) par rapport a l’horizon- 
tale. 

Toute verticale coupant les deux courbes (—C,4) et (—C,) (la premiére 
inversée par rapport 4 son nouvel axe vertical, la seconde retournée par rapport 
a l’axe horizontal demeuré en place) y marque deux points A et B, et deux 


seulement, qui déterminent deux vecteurs A et B dont les composantes hori- 


Vol. VIII, 1957 Comment agit le cable tracteur d’un téléphérique 109 


zontales satisfont a l’égalité Hy + H,; =H tandis que le vecteur vertical BA 
mesurera en grandeur et signe la somme W 4 + Vi = B A. 

Pour qu’une verticale puisse avoir deux points d’intersection A et B avec les 
courbes (—C4) et (—C;,), il suffit qu’elle soit tracée 4 n’importe quelle distance 
de yz supérieure a H si H est positive, et A n’importe quelle distance positive de 


Figure 5a 


Vp Si H est négative. Et les formes de ces deux courbes sont telles qu’une verti- 
cale qui les coupe toutes deux ne peut couper chacune qu’en un seul point. 
Dans ce cas a chacun des points d'intersection correspond un vecteur de traction 


oe 4 pour l’un et un vecteur de traction G; 3 pour l’autre dont les composantes 


horizontales H, et Hg satisfont bien a la relation H, at Fig = —H et dont les 
composantes verticales W, + Vz ont pour résultante géométrique le vecteur 


B Aen grandeur et signe. 

Or si on appelle Z le point commun aux deux courbes (—C,) et (—C,), on 
constate que ce point ne peut pas ne pas exister. Cette existence étant une 
conséquence de l’inégalité des angles 7, et #4 des directions asymptotiques 


(avec #4, > 4). Les verticales bisécantes tracées 4 gauche de Z donnent des 


résultantes BA positives (dirigées vers le bas) dont les valeurs diminuent 
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progressivement de l’infini 4 zéro quand la verticale bisécante se rapproche de 
Z. Quand la verticale bisécante se déplace au- dela de Z, la résultante BA des 
deux composantes verticales des tractions © 4 et e p devient négative (dirigée 
vers la haut) et croit indéfiniment, a partir de la valeur zero, cette croissance 
indéfinie étant la conséquence de l’inégalité des inclinaisons 7, et 7, (avec 7, << 
des deux asymptotes. 


H>0O 


Figure 5b 


Ainsi se trouve démontrée la possibilité de réaliser avec deux brins trac- 
teurs et au point C de leur jonction, n’importe quelle force pourvu que ce point 
C soit compris entre les verticales des appuis A et B du cable tracteur et a une 
. hauteur intermédiaire entre ces deux appuis. 


2° Cas d'un point C a une altitude inférieure & celle de l’appui le plus bas 


Dans ce cas le point C marque Pextremité inférieure des deux sustentes 


fictives AC et CB. Donc les deux tractions C, 4 et C; zp dont la résultante sera 
l’effort exercé par le cable tracteur sur le point C, ont des courbes représenta- 
tives qui traversent l’axe des abscisses avec des asymptotes inclinées (é, et 4) 
de part et d’autre de cet axe. En recommengant le méme raisonnement que 
dans le cas précédent on arrive a la méme conclusion. 
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suspensions réalisant n’importe quelle force < ayant une composante horizontale 
H positive, pour un point C a une altitude intermédiaire entre celles de A et de 


B; la figure 5b montre la méme représentation pour une composante hori- 
zontale H négative. 


La figure 5a montre la représentation de la détermination graphique de deux 


Figure 6a 


Les figures 6a et 6b montrent les mémes représentations pour un point C a 
une altitude plus basse que celle de l’altitude inférieure A. 

Dans les quatre cas types aussi figurés, les composantes verticales V positives 
(c’est a dire dirigées ver le bas) correspondent toujours aux suspensions définies 
par les verticales 4 gauche du point Z, et les composantes verticales négatives 
(dirigées vers le haut) correspondent aux suspensions définies par des verticales 
a droite du point Z. 

Le fait que ces Ss quatre représentations types aient été établies pour des 
sustentes fictives AC et CB ayant des projections horizontales égales (verticale 
du point C au milieu de l’intervalle AB) n’infirme nullement la conclusion du 
cas général puisque la représentation de l’une ses deux courbes (—C4) ou 
(—C,) y deviendrait modifiée homothétiquement a partir de son foyer y4 ou yz. 

Ainsi seule l’intervention de deux brins tracteurs permet de réaliser en tout 
point C accessible au chariot du téléphérique telle force que l’on voudra sans 
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restriction aucune, hormis celle qu’imposera la limite de résistance du cable 
tracteur. Et il ressort péremptoirement de cette étude que dans leur action 
conjuguée de traction aucun des deux brins n’assume un réle prépondérant sur 
celui de l’autre: ils sont rigoureusement indispensables l’un a l’autre, et la 
défaillance de l’un entraine automatiquement l’inutilité de l'autre et ’impotence 
du téléphérique. 


Figure 6b 


x 


5. Restriction dans le choix de la force 4 réaliser en C du fait de la 
limite de résistance 4 la rupture des brins du cable tracteur 


L’étude de cette restriction se fait trés simplement en utilisant la courbe 
inversée (—C,4) par rapport a la verticale et la courbe retournée (—C,) par 
rapport a l’axe des abscisses, mais en tracant ces courbes de telle sorte que leur 
asymptote verticale soit commune. On SOS ainsi une représentation du 
genre de celle de la figure 7. 

A partir du foyer y commun de ces deux courbes, nous tracerons une cir- 
conférence de rayon Ly, dont le rayon mesure en longueur la résistance a la 
rupture (ou une fraction de sécurité de celle-ci) du brin issu de l’appui A. Cette 
circonférence coupera la courbe inversée (—C,4) en deux points A’ et A” sis de 
part et d’autre du point A,, s'il y a eu intersection: un allongement suffisant 
de la sustente fictive AC, dont l’inclinaison 7, demeurerait inchangée, pourrait 
rendre A, égal ou supérieur a Ly. 

On opérera de méme avec un autre rayon Lz pour tracer du foyer commun y 
une autre circonférence qui coupera la courbe retournée (—Cg) en deux points 
B’ et B” sis également du part et d’autre du point B,,. Le rayon Lp sera rarement 
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égal 4 Ly pour deux raisons: les deux brins ne présentent pas nécessairement 

la méme résistance a la rupture et méme s’il en était ainsi il faudrait retrancher 

a celle du brin B le poids d’une longueur de cable égale a la dénivellation entre 

Vappui B et le point C étudié. 

. Si l’on trace ensuite les 4 verticales passant par les 4 points susdéfinis: 
A’, A”, B’, B”, on voit que le choix de la force A réaliser sera limité par un 


Figure 7 


éloignement l’une de l’autre des deux asymptotes (des deux foyers y4 et yz) 
qui fera coincider la verticale de B” avec celle de A’ et par un autre éloignement 
amenant en coincidence la verticale de B’ avec celle de A”. Dans le second cas 
la composante horizontale H de la force limite aura une valeur maximale 
H}, + H", (H’, valeur positive maximale H) valeur négative, minimale en valeur 
absolue) tandis que sa valeur minimale Hj, + H’ sera atteinte dans l'autre cas 
extréme d’éloignement. Dans l’un comme. dans l’autre cas les deux brins 
tracteurs seraient sollicités 4 leurs limites respectives de résistance. 

Si de plus on considére tous les cas admissibles d’éloignement entre les deux 
positions extrémes résultant des limites respectives de la résistance a la rupture 
des deux brins tracteurs, si on considére encore que des modifications de lon- 
gueur des sustentes fictives a inclinaisons inchangées, interviennent également, 


_ ZAMP VIII/8 
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on en déduira qu’il est impossible d’affirmer a priori qu’un brin tracteur sera 
systématiquement moins sollicité que l'autre. 

Si on pourvoit un des brins tracteurs d’un dispositif lui imposant une 
traction fixe en un point donné (qui sera pratiquement le point fixe d’ot il est 
issu) il n’en résulte nullement une garantie de limitation pour la traction dans 
l’autre brin; seul sera restreint le choix de la force pouvant étre réalisée par la 
conjonction en C des deux brins. 

En effet supposons que cette traction fixe soit imposée au brin issu de B 
et que celle-ci soit figurée sur la courbe représentative des tractions Cz par 
deux points B’ et B” (sur une circonférence de centre y,) et de coordonnées 
respectives Hy, Vz et H%, Vz. Dans le choix de la force réalisable, la compo- 
sante H de celle-ci devra déja étre choisie telle qu’aprés avoir déplacé l’asymp- 
tote verticale de la courbe (—C,4) de la quantité H l’une au moins des deux 
verticales tracées par les points B’ et B” se trouve a droite de l’asymptote 
verticale déplacée y,. On en déduit que H ne peut avoir qu’une valeur égale ou 
inférieure a celle de Hz (puisque Hz > Hz). 

Pour les valeurs de H comprises entre Hz et Hz on ne dispose que d’une 
seule verticale bisécante, celle qui passe par B’, donc d’une seule valeur de V 
par valeur de H; ce ne sera que pour les valeurs de H inférieures 4 HZ que l'on 
disposera de deux verticales bisécantes, celles passant par B’ et B”, donc de deux 
valeurs de V par valeur de H. 

Il suffit, a l'aide de la figure, d’imaginer par la pensée l’infinité de combi- 
naisons possibles (inclinaisons des sustentes de leurs longueurs relatives) pour 
conclure qu’un dispositif imposant une traction fixe a l’un des brins ne donne 
aucune garantie quant a la limite de l’effort qui pourra étre imposé a l’autre brin. 

Si cet autre brin est aussi pourvu d’un dispositif analogue lui imposant une 
traction fixe, il devient évident que en chaque point C le choix de la force 
réalisable par la conjonction des deux brins tracteurs deviendra limité a 
4 vecteurs forces, ceux résultant de la mise en alignement vertical des deux 
points limites A’ et A” avec les deux autres points limites B’ et B”. 

On a coutume de qualifier de dispositif de réglage de tension celui d’un contre- 
poids transmettant une traction constante au brin qui contourne une poulie. Cette 
dénomination est inexacte; il ne s’agit pas d’un réglage mais bel et bien d’un 
dispositif imposant une tension, en un mot d’un dispositif de mise en tension pour _ 
employer la dénomination dont nous nous sommes servis au début de cette étude 
quand nous avons donnée la description générale schématique d’un téléphérique. 

Le véritable dispositif de réglage de tension est celui réalisé par la molette 
motrice, qui regoit du moteur des impulsions tractrices vraiment réglables. 
Aussi du point de vue de la sécurité des cables tracteurs la meilleure organi- 
sation d’un téléphérique est celle qui fait du brin rétroactif le brin de mise en 
tension (grace a l’action d’un contrepoids) et du brin actif le brin de réglage 
de la tension (grace a l’action modulée du moteur). 
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En limitant a cet exposé le développement de quelques déductions parmi les 
trés nombreuses que permet la confrontation des courbes représentatives (A) et 
(B) des tractions dans un cable a ses extrémités nous pensons avoir présenté la 
solution de ce probléme annexe que nous avions invoquée dans un mémoire 
antérieur pour proclamer l’impropriété de la dénomination cAble-lest. Certes 
notre affirmation a été contestée, curieusement sans aucune démonstration, 
car pour certain la critique de la conclusion d’un probléme est plus aisée que 
art d’en trouver la solution. 

Nous ferons remarquer enfin que pour de telles études l’analyse ne peut 
fournir des éclaircissements comparables 4 ceux qu’offre le raisonnement, non 
pas que les mises en équation soient difficiles, mais seulement parce que celles-ci 
aboutissent a des relations entre fonctions et arguments hyperboliques dont 
l’interprétation devient vite inextricable. 


Annexe 


Détermination de la forme et des caractévistiques structurales des courbes (A) et (B), 


veprésentatives des tractions A et B exercés par un cable de longueur variable fixé 
en deux points A et B 


On sait que toute caténe formée par un cable sous une sustente AB est 
obtenue par une transformation homothétique d’un certain arc de chainette 


_—. 


A,, B;, qu'une corde A,, B,, ayant méme inclinaison / que la sustente, découpe dans 
une chainette de référence, le rapport d’homothétie étant 4 B/A, By. 

Prenons comme chainette de référence celle dont la corde en position minimale 
d’inclinaison i a pour longueur celle de la sustente 4B. Le module 7 de cette 
chainette, ainsi que les ordonnées canoniques correspondantes A a et Bb de la 
sustente AB sont déterminées par les caractéristiques de la position minimale en 
inclinaison 7. 

Ainsi chaque caténe de suspension sous la sustente A B est engendrée par une 
certaine corde A, B,, de cette chainette de référence, Mais comme les cordes 
A, By peuvent ves appareillées deux a deux, l’une en dessous de la corde mini- 


ale AB (et nous l’appellerons Aj Bj) l'autre en dessus (et nous l’appellerons 
A‘, By), de telle sorte que A, ai/A’, BY == Aah was , BY, il existe toujours deux 
suspensions, l’une lache l'autre tendue, piiduiaet & chacune des extrémités de 
Ja sustente des tractions Ac égales mais de directions différentes, de méme que les 


tractions By: 
Nous appellerons 7; le module de la suspension tendue et 7), celui de la sus- 


pension lache. 
Nous appellerons V;, H;, les composantes, verticale et horizontale, de la trac- 
tion A, en suspension tendue, et V/ et Wj ces mémes composantes en suspension 
it 1 k k 


lache. 

Pour les composantes des tractions B, a lextremute superieure B, nous em- 
ploierons les lettres W, et W, pour distinguer les composantes verticales en 
suspension tendue et lache, i composantes horizontales demeurant respec- 
tivement égales a celles des tractions AL au signe prés puisqu’elles sont toujours 
orientées vers la zone 4 l’aplomb sous la sustente. 
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Convenons de donner une valeur positive aux composantes verticales V et W 
quand elles sont orientées versle bas (comme la pesanteur) et convenons de 
compter H positivement quand cette composante horizontale est dirigée de 
l’aplomb de l’appui le plus bas vers celui de l’aplomb le plus haut. (Au cas ou les 
appuis seraient a la méme altitude, l'étude est si simplifiée qu’une représentation 
graphique des courbes (A) et (B) devient superflue.) 


Ces notations étant ainsi précisées, on en déduit les relations suivantes: 


vy! v A B W Y A B " 
ko k 7 z k Sa a k? 
AUB Are 
A’a,  __. Ajay, 
[igs AN Opn awn! Me) AB —~=AB4q,, 
/, "2 Ae yet 
kk k k 
“ B’ bf B!" b! 
B, = Bb, = AB + = AB(q, + sini), 
7 Ee AUR 
ES ae hess es 
ee la 7 AB, pe a le 
A, B, Aj, Bi, 
Re (S 
Wi = AB, Whe AB 
ALE At B" 


Examinons comment évoluent les suspensions engendrées par des cordes 


génératrices qui s’éléveraient progressivement dans la chainette de référence, a 
partir du point J de celle-ci ot la tangente a V’inclinaison 7. 


t 
‘ 
y 
1 


: ae | 
‘ ‘ 
: _ 


a Gap OL Xa Me Selb. 10h 


Figure 8 
%, point de rencontre des tangentes aux extrémités de la corde AB en position minimale pour 


Vinclinaison i (A a | ABaune valeur minimale ainsi que B b / AB); %y, point de rencontre des aligne- 
ments des extrémités de deux cordes paralléles en correspondance de traction pour V’inclinaison 7. 


eee Apa, et Bh bk _ _ Bho 
An Bh Aj, BE Ap Be Ah, By 


Vol. VIII, 1957 Comment agit le cable tracteur d’un téléphérique Ny 


En suspension infiniment tendue (corde génératrice infiniment voisine du 
point J) les deux tractions (A) et (B) sont infiniment grandes, égales et apposées 
suivant la direction de la sustente. 

Au fur et 4 mesure que la suspension devient moins tendue, les tractions A 
et B commencent par devenir de plus-en-plus petites, A décroissant plus vite que 
B (le rapport B/A augmente). 

Cette décroissance simultanée de A et de B se poursuivra jusqu’a ce que la 
corde génératrice A, B, se confonde avec la corde minimale 4 B; puis cette corde 
génératrice porunteant ds’élever, A et B se mettront a croitre, A croissant plus vite 
que B (le rapport B/A diminue et tend vers 1). Les composantes horizontales (H CU 
—H) qui étaient infinies ne cessent de décroitre et cette décroissance se poursuivra 
jusqu’a la valeur nulle, qui correspond 4 une suspension infiniment lache en- 
gendrée par une corde génératrice 4, B,, tracée a l’infini au-dessus du point J. 

On notera que pendant cette Svoliitienk langle B que fait la diel aet de la 
traction B avec l’horizontale n’a cessé de croitre pour tendre vers 7/2, la traction 
étant toujours orientée vers le bas. 

Au contraire pour la traction A, d’abord orientée vers le haut, on voit que sa 
direction s’abaisse progressivement jusqu’a devenir horizontale (la corde généra- 
trice est alors celle qui passe par le sommet de la chainette de référence), puis cet 
angle « se remet a croitre pour tendre vers 7/2. Il suffit d’attribuer a4 cet angle « 
une valeur négative pour les tractions A orientées au-dessus de l’horizontale, pour 
pouvoir conclure que cet angle a, ainsi mesuré algébriquement, ne cesse de croitre 
pour tendre vers 2/2 comme l’angle . 

Il est évident que pour une caténe quelconque de suspension la composante 
verticale W d’une traction B, a l’extrémité la plus haute, est toujours supérieure, 
en valeur absolue, a la composante verticale V de la traction a l'autre extrémité. 

Si l’on compare les 4 composantes verticales de deux suspensions appareillées, 
on constate que tant que les cordes génératrices A, Le i ee BY sont voisines de 
la corde minimale on a W, > V;, mais qu’il existe fe Spivciiane et un seul 
poo lequel W, = Vj, et qu’au ree de cet appareillage W, est toujours inférieur 

a V;. 

ia relation qui donne la valeur V; en fonction de la corde génératrice Aj, Bi, 

de la suspension correspondante, soit: 


A; S ra 
Ay B; 


Vv, = 


nous permet de démontrer: 
lo que la valeur absolue de V, ne cesse de décroitre quand on passe de la 


suspension infiniment tendue (pour laquelle elle a une valeur infinie), a la sus- 
pension — encore tendue — a tangente horizontale a l’extrémité inférieure (corde 
génératrice S B,) auquel cas V, a une valeur nulle; 

20° qu’a partir de cette suspension tendue a tangente horizontale, la compo- 
sante V, prend des valeurs positives croissant indéfiniment au fur et 4 mesure que 
la suspension devient moins tendue, puis minimale, puis lache. 

Pour cela nous aurons recours a un raisonnement géométrique différentiel 
appliqué aux variations des longueurs de deux cordes génératrices infiniment 
voisines et de celles de leurs arcs. 

Considérons d’abord une corde génératrice asuspension plongeante (c’est-a-dire 


une corde qui coupe l’axe vertical de la chainette). Soit A(Ax eS) le petit élément 
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d’arc qui sépare l’extrémité inférieure de la corde A,, B; de celle infiniment 
voisine. Ce petit élément d’arc peut étre confondu avec |’élément rectiligne de 


méme longueur aligné suivant la tangente A; 7;,. Si la nouvelle corde infiniment 
voisine est au-dessus de A, B,, sa longueur aura été accrue d’une quantité 
A(A, B,); Vélément A(A;, S) et celui A(A;, B;,) sont les deux cétés d’un triangle 
semblable au triangle A, T;, B,, donc on a la relation 


A(4y S), An Te 
A(A;, B;) Ay B,, 


Mais le point de rencontre T;, des deux tangentes a la chainette aux extrémités 
de la corde A, B, (point T,, qu’il paraitrait opportun d’appeler «péle de la corde 
A, B,) est dans la moitié verticale de la chainette opposée a celle du point A,. 
Donc la développée a cette chainette, a partir du sommet et vers le coté ot se 


trouve A,, coupe le vecteur A, 7, en un point compris entre A, et T,, donc 


SS 


A, TI; est plus grand que l’are A, S. 
On en déduit que A(A, S)/A(A, B,) est supérieur a A, S/A;, By. 
Pour savoir ce qu’est devenu V, au cours d’un déplacement infinitésimal vers 
le haut de la corde génératrice A;, 6, on comparera 


Figure 9 


A’ A" = B’B" = A(AB); AA’ = AAS); BB’ = A(BS). 
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Or quand on augmente les deux termes » et @d’une fraction njd de quantités 
respectivement égales a n’ et A d’, dont le rapport n’/d’ est supérieur A n/d, la 
nouvelle fraction augmente de valeur. 

Done V,, augmente de valeur quand sa corde génératrice s’éléve; autrement dit 
quand le module de la suspension diminue, 

Pour les cordes génératrices dont l’extrémité inférieure est du méme cété que 
B,, par rapport a l’axe, la croissance de la valeur absolue de V au fur et A mesure 
que la suspension devient de plus en plus tendue est évidente — sans qu'il soit 


. . ce . . . , . . eee 
besoin de recourir 4 un raisonnement géométrique différentiel — puisque A, S 


augmente (et tend vers I S) tandis que A , B, décroit. 
L’étude des composantes verticales 


procéde du méme raisonnement géométrique différentiel, lequel fait conclure que 
W, toujours positive, a une valeur minimale pour la suspension engendrée par la 
corde 4, B;, dont le péle J; (point de rencontre des tangentes a la chainette aux 
deux extrémités) se trouve sur la développée a la chainette issue du sommet de 
celle-ci et tracée vers le coté ot se trouve le point B,. 

Ce point /, existe toujours car on démontre que la polaire caténique d’une 
direction 7 (lieu des points de rencontre des tangentes aux extrémités des cordes 
paralléles d’inclinaison 7) est une courbe qui ne peut pas avoir de tangente hori- 
zontale, ni aucun point d’intersection avec aucune autre de méme nature. On 
en déduit notamment que quand le point J de la chainette ot la tangente a 
Vinclinaison @ (lequel point 7 marque donc l’extrémité supérieure de la polaire 
caténique de la direction 7) se déplace sur la chainette, le point J de sa polaire 
se déplace dans le méme sens sur la développée de la chainette. 

Parmi ces points J il convient de distinguer celui unique, Jy, ot la tangente 
a la développée est aussi tangente a la chainette. Les coordonnées de cette corde 
remarquable se calculent aisément a partir des deux équations 


ree a! ghee Sih shi tes 2k 
Y y y Y 
dot ; 
cat a 7s ek —argsh {/ ——— 
- sh “1 
e 


dont la racine est *,/y = 1,20938 d’ot l’on déduit avec une précision de 1/10000 
%/¥ = — 0,61548. 

Avec une précision de 1/200 l’extrémité supérieure de cette corde remarquable 
peut étre confondue avec le point de la chainette ou la tangente passe par le point 
de rencontre de l’axe et de la base 


ee eh LOVE 1200.38 
V4 


Il suffit de considérer la figure ot est représentée la tangente commune EEG Mi; 
_ pour se rendre compte que toutes les autres cordes ayant leur pole caténique sur la 
développée auront leur extrémité inférieure kK, comprise enire Aowiet pon ete rd: 
chaque point K, correspondent deux pdles caténiques J, et J, appartenant a des 


: % : ar A Al 4 
_ polaires issues de point J, et I; (avec 7% < %< 7{). On remarquera accessoirement 
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que l’angle en J, est obtus tandis que celui en J,’ est aigu. Ainsi se trouve démontré 
que la composante verticale minimale de la traction a l’extrémité supérieure 
d’une caténe de suspension est obtenue pour une caténe tendue, comprise entre 
la caténe minimale et celle 4 tangente inférieure horizontale. 


TMs 
Figure 10 


Un autre élément important de la structure des courbes (A) et (B) est leur 
asymptote, figurative des tractions aux extrémités des caténes de suspension 


infiniment tendues. 

Ces asymptotes font avec l’axe des abscisses H un angle égal a celui, 7, que 
fait la sustente AB avec l’horizontale. La mise en place de ces asymptotes sera 
achevée quand nous aurons démontré que chacune coupe l’axe des ordonnées a 
une distance de l’origine telle que ce vecteur mesure le poids d’une longueur de 
cable égale 4 A B/2. 

En effet, en comptant positivement les arcs de la chainette de référence, dans 
le sens de gauche a droite (l’extrémité inférieure A étant supposée a gauche), on | 


a en grandeur et signe: 


Ap A pa aes 
Vi= + ALS = + (4;1+75), 
A ee Be 
WwW, = BL S= = Z (Bil +28) 
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La relation qui définit la composante horizontale H,, d’une traction A, ou B,. 
(composante horizontale mesurée positivement-de gauche a droite) est 


Hy y 
Ai Bi 


pour toutes les tractions B,. de l’extrémité supérieure et 


pour toutes les tractions A. de l’extrémité inférieure de la sustente. 
Tenant compte de ces deux relations, on peut donner aux expressions de V,, 
et de W, les formes suivantes: 


ou H; est la composante horizontale de A, (donc toujours positive) dans l’équation 
de V; et la composante horizontale de B, (donc toujours négative) dans 1’équation 
de W,. 

A la limite, quand la corde génératrice A; Bj, se rapproche de plus en plus du 
point J, ot la tangente a l’inclinaison 7 de la sustente, Arias B,, tend vers I/2, 
tandis que B, T/A; B/ tend — 1/2. 

D’autre part, I S/r est égal A — tgi. 

On a donc a la limite 


Va HY tet et We Gaon ee ni tei: 


Et par des raisonnements analogues, ne procédant que de considérations géo- 
métriques, on démontre aussi que les deux courbes (A) et (B) ne peuvent couper 
leur asymptote, ce qui achéve de donner de la structure de ces courbes toutes 
les précisions suffisantes pour fonder sur elles l’étude que nous venons d’exposer 
concernant les cables tracteurs. 


Summary 


Many studies, and even official regulations use the name ‘balance cable’ to 
describe one of the two ends of the traction cable of a telpher. 

This definition implies that the end of the cable referred to as balance cable 
plays a minor part in the efficient running of the telpher. 

The author shows this interpretation is not justified, in the sense that the 
two ends of the cable are mechanically equivalent and indispensable to the 
system. The article shows also that the counter-balance cannot be used as a 


device adjusting stress during operation, this effect being obtained only throug 


the action of the driving pulley. 


(Recu: le 11 aotit 1956.) 
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Géométrie de la chainette ou hypergéometrie 


Par Marie Marcet Nicoras, Paris?) 


La solution d’un probléme peut étre obtenue soit par la méthode analytique, 
soit par le raisonnement ; la premiére montre la solution, la seconde la démontre. 
Une solution n’est parfaite que si elle a été démontrée. 

La méthode analytique appliquée a l'étude des systémes comprenant des 
cables est généralement stérile parce qu’elle fait intervenir des relations entre 
des fonctions hyperboliques et leurs arguments et parce que rares sont celles de 
ces relations qui soient susceptibles d’une interprétation facile. 

Le raisonnement demeure malaisé tant que l’on ne dispose pas d’une 
géométrie qui, par son essence, convienne a ce genre de problémes. 

En effet tout probléme de mécanique ot interviennent des cables comporte 

toujours la confrontation de longueurs de cables suspendus en des points 
donnés de l’espace avec les tractions qu’il convient d’exercer sur eux en ces 
‘points de fagon a satisfaire 4 certaines conditions imposées. Ces cables, sup- 
posés infiniment souples, mais non dépourvus de masse, se trouvent soumis a 
un champ de gravitation que l’on peut admettre d’intensité constante et de 
direction uniforme dans le cas des problémes usuels. 

Sous l’influence d’un tel champ, celui de la pesanteur, les cables prennent 
la forme d’arcs de chainette, de caténes, découpés sur des chainettes dont l’axe 
est paralléle a la direction du champ, donc a la verticale. 

La longueur des caténes ainsi tracées entre deux points A-et B de l’espace 
est fonction non seulement de la distance mesurée rectilignement entre ces 


points, mais aussi de l’inclinaison du vecteur AB au-dessus du plan horizontal 
et enfin de la traction exercée sur le cable en un de ces points, traction qui 
détermine le module de la chainette sur laquelle est découpée la caténe, et 
vice versa. 

On congoit combien est inadaptée a l'étude de tels problémes une géométrie 
fondée sur l’emploi de la régle et du compas, c’est-a-dire sur une mesure rigide 
et rectiligne des distances. 

Dans un monde soumis a une gravitation d’intensité infinie — ou ce qui 
revient au méme dans un monde dont tous les matériaux seraient infiniment 
souples — les besoins de mesure et les moyens matériels de les satisfaire auraient 
imposé a des étres raisonnables la conception d’une autre géométrie que celle 
usuelle. Cependant si les besoins et les moyens de mesure des étres raisonnables 


1) Secrétariat d’Etat aux Travaux Publics, aux Transports et au Tourisme. 
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vivant dans cette ambiance physique demeuraient limités’ 4 un univers 
perceptible infiniment petit, c’est-a-dire dont la plus grande dimension ne 
dépasserait pas le millioniéme ou le milliardiéme du plus petit module réalisable 
sur les régles caténaires de mesure employées par ces étres, les apparences 
leur feraient encore concevoir une géométrie euclidienne identique a celle dont 
nous nous servons. 

Mais ceux-ci commenceraient a douter de l’universalité de cette géométrie 
au fur et 4 mesure que s’étendrait la portée de leurs perceptions. 

Le méme doute devait s’emparer des géométres usant des régles rigides 
quand ils substituérent a celles-ci des fils constitués par des rais lumineux 
soumis a des champs de gravitation. 

Or l'étude des problémes de cables, méme envisagés a la trés petite échelle 
des besoins de la technique, impose |’emploi de cette géométrie ot les distances 
sont mesurées suivant des caténes. 

Nous nous proposons d’exposer les fondements de cette géométrie speciale. 
Comme dans le cadre de ces Annales il serait impossible d’en faire un exposé 
complet, nous nous bornerons a décrire et a étudier sommairement celles des 
lignes essentielles au recours desquelles nous ont conduit nos études sur les 
cables de téléphériques. 

Pour la commodité de l’exposé nous appellerons cette géomeétrie «hyper- 
géométrie» rappellant ainsi qu’elle traduit des relations entre des fonctions 
hyperboliques. Pour faciliter la confrontation de la géométrie usuelle avec’ 
Vhypergéométrie nous attribuerons aux lignes de cette derniére, qui tiennent 
le méme -réle qu’en géométrie euclidienne, les mémes désignations respectives 
précédées du préfixe «hyper»; nous aurons donc des hyperdroites, des ee 
circonférences, des hyperspheres, des byprperaileies, des hypernormales, etc.. 
etc... et méme des hyperhyperboles. 

Cette hypergéométrie en suggére d’autres, adaptées a des champs de gravi- 
tation qui n’auraient plus des intensités et des directions (au sens euclidien) 
uniformes mais qui correspondraient par exemple a des champs de gravitation 
concentriques, au sens usuel comme au sens hypergéométrique du terme. 


Droite et hyperdroite 


La droite est la ligne selon laquelle on mesure la distance entre deux de ses 
points. L’hyperdroite sera donc une chainette d’axe vertical, plus exactement 
d’axe paralléle 4 la direction uniforme du champ de gravitation d’intensité 
constante auquel sont soumises les régles de mesure supposées infiniment 
souples. L’hyperdroite est donc toujours contenue dans le plan vertical défini 
par les deux points dont on veut mesurer l’hyperdistance. 

Le module de la chainette de l’hyperdroite est fonction de la tension de la 
régle souple servant 4 mesurer les hyperdistances. L’hyperdistance entre deux 
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points dépend donc de la distance de ceux-ci, de l’inclinaison de leur alignement 
au-dessus de l’horizontale et de la tension de la régle. 

La différence fondamentale entre la droite et l’hyperdroite est que les 
points de la premiére n’ont pas de personnalité géométrique tandis que ceux 
de la seconde en ont une. 

Par exemple la direction de la tangente en un point de l’hyperdroite est 
fonction du point choisi sur celle-ci; de cette remarque il sera tiré parti dans 
l’étude de Vhyperplan. 


Figure 1 


e468 AB AG 4 CB SABA ACERS Ae. 


Mais il existe d’autres dissemblances remarquables entre la droite et 
Vhyperdroite. 

En effet alors qu’entre deux points on ne peut tracer qu’une seule droite, 
on peut toujours tracer deux hyperdroites équimodulées entre deux points 
(l'une est une caténe soumise a des efforts internes de traction, l’autre est la 


méme caténe symétrique de la précédente par rapport au milieu de AB et 
soumise a des efforts internes de compression) (figure 1). 

La ligne droite est la ligne suivant laquelle la distance mesurée entre les 
deux points A et B est la plus courte; plus exactement tout point C pris en 
dehors de la droite AB détermine deux éléments de droite AC et CB dont la 
somme des longueuts est supérieure a AB. 

Les deux hyperdroites équimodulées passant par A et B délimitent deux 
zones dans le plan. Tout point C, de la zone externe détermine deux hyper- 
droites équimodulées (dans le module considéré) telles que la somme des 


hyperdistances AC, + Cab. soit supérieure a V’hyperdistance AB. Tout point C; 
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de la . zone interne détermine deux x_segments dhyperdroites équimodulées tels 


que AG + C.B B soit inférieure 4 AB mesurée suivant l’hyperdroite. 

On en déduit que la propriété fondamentale des triangles, savoir que la 
longueur d’un cété est inférieure 4 la somme des longueurs des deux autres 
cétés, n’est plus une propriété absolue des hypercétés d’un hypertriangle. Il 
suffit en effet qu’un sommet se trouve 4 l’intérieur de la zone délimitée par un 
hypercoté et le symétrique inversé de celui-ci pour que la longueur de cet 
hypercoté soit supérieure a la somme des deux autres; elle lui serait égale si le 
sommet était sur le symétrique inversé. 


Plan et hyperplan 


I] existe du plan deux définitions: Le plan est la surface qui contient toute 
ligne droite joignant deux points quelconques de cette surface. Le plan est 
aussi la surface engendrée par la rotation d’une ligne droite autour d’un axe 
qui lui est orthogonal. 

La premiére définition nous fera donc appeler hyperplan tout plan vertical, 
car seule une telle surface contient toutes les hyperdroites joignant deux 
quelconques de ses points. 

La seconde définition nous ménera a la méme conclusion quand nous 
aurons éliminé le cas d’une rotation autour d’un axe orthogonal a une hyper- 
droite au sens restreint, au sens ponctuel, du mot. Il faut en effet que, comme 
dans le cas de la droite, l’axe de rotation soit orthogonal a tous les éléments 
rectilignes infiniment petits qui constituent l’hyperdroite et non pas a un seul 
d’entre eux. 

Or la seule direction orthogonale a tous les éléments rectilignes infiniment 
petits d’une hyperdroite est celle qui est normale au plan vertical qui contient 
cette chainette. 


Circonférence et hypercirconférence 


La circonférence est le lieu des points d’un plan qui sont a une distance 
donnée d’un point de ce plan. Les lignes droites qui servent 4 mesurer cette 
distance peuvent étre considérées comme les arcs d’une chainette infiniment 
tendue, donc comme des caténes ayant méme module, savoir un module de 
valeur infinie. 

Tenant compte de cette derniére remarque nous définissons hypercircon- 
férence comme étant le lieu des points d’un hyperplan (plan vertical) qui sont a 
une distance donnée d’un point de cet hyperplan, étant entendu que cette 
distance est mesurée le long d’arcs découpés sur une chainette d’un module 
donné (et d’axe vertical évidemment). Cette distance et ce module définissent 


Vhyperrayon. 
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On voit immédiatement que les hypercirconférences concentriques ayant 
toutes méme longueur d’hyperrayon (mais tracées sous des modules variables) 
sont toutes enchdssées les unes dans les autres, celles de modules supérieurs 
entourant celles de modules inférieurs. Les formes extrémes de ces hyper- 
circonférences, ayant méme longueur d’hyperrayon, sont d’une part la circon- 
férence de rayon J, d’autre part le vecteur vertical de longueur 2/ (figure 2). 


_-0,0921 X 1 (i/m = 100) 
~~ 0,3 X 1 (I/m = 20) 
_-—0,881 X 1 (l/m = 1) 


Figure 2 


Hypercirconférences d’hyperrayon /, constant mais tracées sous divers modules. 


ee [ec ie 
m 


OH 


1 


1 0,98964 | 0,96248 | 0,88137 | 0,7966 | 0,7217 0,3 0,0921 0 


Pour que l’hypercirconférence se distingue pratiquement de la circonférence, il faut que son hyper- 
rayon soit supérieur au module sous lequel elle est tracée. 


Toutes les hypercirconférences concentriques d’hyperrayon / ont en commun 
les deux points situés sur la verticale du centre a une distance / de celui-ci. 

Nous appellerons diamétre vertical d’une hypercirconférence le vecteur 
joignant les deux points de celle-ci situés sur la verticale passant par le centre, 
et diamétre horizontal le vecteur horizontal passant par le centre et limité par 
Vhypercirconférence. 

Rapportée a son centre et 4 des coordonnées cartésiennes, u et v, dont les 
axes sont suivant le diamétre vertical et le diamétre horizontal, l’hypercircon- 
férence d’hyperrayon J, mesuré sous le module m, a pour équation: 


v+4m% sh?“ =]2, 
2m 


En effet entre les abscisses canoniques % et x, de deux points d’une chai- 
nette de module m et délimitant sur celle-ci un arc de longueur J, on a les 
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relations Z 
m (sh + — sh ==.) =(=—2m sl (= ed (a= * 
ra rh sh ae ch oo 
= 2msh | (* i 
om \ on 
e€ 
m (ch —ch =o) =v=2msh—"— sh (+ +o) 
es m m 2m ; 2m 
ou 
72 ek 
4m? sh? 4m? tie. i 
2m 2m 
et 


¢ Uu 
J? — y? = 4. m? sh? 
2m 


L’étude d’un arc de longueur constante tracé sur une chainette montre que 
la projection horizontale de cet arc va en augmentant en méme temps que 
diminue la longueur de sa corde quand diminue Vangle que fait sa corde avec 
Vhorizontale. 

Done le rayon vecteur OP d’un point P d’une hypercirconférence va en 
augmentant quand augmente l’angle que fait ce rayon vecteur avec l’hori- 
zontale. La valeur minimale du rayon vecteur est celle du rayon horizontal, 
égal a: , 

2m arg ere ’ 


La propriété fondamentale de l’hypercirconférence concerne la direction de 
sa tangente en un point. En tout point d’une hypercirconférence la tangente 
est paralléle a la bissectrice des deux tangentes aux extrémités de lhyperrayon 
du point. 

En effet soient AT et BT les deux tangentes d’une chainette aux points 
A et B délimitant un arc de longueur /. Un arc de méme longueur, avec extré- 
mités infiniment voisines de celles du premier, aura ses extrémités A’ et B’, sur 
les deux tangentes AT et BT; si A’ est compris entre A et T, l'autre nouvelle 
extrémité, B’, sera au-dela du point B sur le vecteur TB (figure 3). 


Si par B’ nous menons un vecteur B’ B” equipollent au vecteur A’A, le point 
B" sera un point de l’hypercirconférence de centre A et d’hyperrayon 


fA Bice A'B?. L’égalité de ces deux arcs infiniment voisins implique que 
A’A = BB’, donc a la limite I’élément BB” de Vhypercirconférence est bien 
orthogonal a la bissectrice de l’angle extérieur 4 ATB, donc est paralléle a la 


bissectrice de cet angle. 
L’étude des arcs de longueur constante tracés sur une chainette montre que 


leur bissectrice tangentielle (bissectrice de l’angle formé par les deux tangentes 
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aux extrémités de l’arc) passe progressivement de la direction verticale a la 
direction horizontale quand la corde de l’are passe de la direction horizontale 
a la direction verticale. 

Cette étude montre aussi que quand la corde de l’arc de longueur constante 
passe progressivement de la position horizontale a la position verticale la 
bissectrice tangentielle, d’abord normale a la corde horizontale, ne cesse de 


SS ee 


Figure 3 
AVAL =) BB US wBOB 


former avec la corde un angle aigu, pour atteindre 4 nouveau la valeur 7/2 
quand la corde est devenue verticale. 

On en déduit que la tangente a l’hypercirconférence s’incline progressive- 
ment de la verticale a l’horizontale quand le rayon vecteur OP pivote de la 
position horizontale a la position verticale. 

Toute hyperdroite passant par le centre d’une hypercirconférence ayant 
le méme module est coupée par celle-ci en deux points dont l’hyperdistance est 
constante; celle-ci mesure la longueur de tous les hyperdiamétres de l’hyper- 
circonférence. 

L’hypercirconférence est susceptible d’une interprétation géométrique. Si 
on trace dans un plan vertical la courbe 


U 


et si on projette cette courbe sur un plan horizontal passant par la ligne des 
abscisses u de telle sorte que le vecteur projetant satisfasse 4 la double con- 
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dition d’avoir une longueur constante / et d’dtre orthogonal a la ligne des 
abscisses “, on obtient l’hypercirconférence (figure 4): 


12 — y2=4m%sh2_ 
2m 


------------------3 


in 
we 
, 


Figure 4 


Si on considére tous les points des hypercirconférences concentriques de 
modules variables mais d’hyperrayon constant, qui sont alignés sur le centre 
commun, les tangentes en chacun de ces points aux diverses hypercirconférences 
ont des directions qui s’inclinent progressivement de la verticale a la direction 
orthogonale a celle de l’alignement quand on passe de l’hypercirconférence de 
module zéro (vecteur vertical) 4 |’hypercirconférence de module infini (circon- 
férence de rayon /). 

En effet cette succession de caténes de longueur constante, a4 modules 
croissants est caractérisée par le fait que le rapport de la longueur de la sustente 
(corde de la caténe) a la longueur de la caténe sustendue va en augmentant en 

tendant vers l’unité. 

Puisque ce rapport va en augmentant cela signifie que dans une chainette 

de référence, de module fixe, la corde de la caténe figurative de l’hyperrayon 


ZAMP VIII/9 
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se déplace parallélement a elle-méme en se rapprochant du point ot elle devient 
tangente a la chainette; la bissectrice de l’angle formé par ses tangentes s’incline 
alors progressivement de la verticale a la direction orthogonale a celle dela corde. 
Cette propriété est d’ailleurs confirmée par l’analyse. 
En effet en écrivant que la différentielle de /? est nulle, on trouve pour 
l’inclinaison dv/du de la tangente 4 une hypercirconférence de module m: 


ou apparait une des rares relations entre une fonction hyperbolique et son 
argument, dont l’interprétation ne présente aucune difficulté. 

Deux hypercirconférences sont homothétiques quand leurs hyperrayons et — 
leurs modules sont dans le méme rapport. 

La condition nécessaire et suffisante pour que deux hypercirconférences | 
soient homothétiques est que la droite qui joint leurs centres les coupe en des 
points ot les tangentes soient paralléles. 

Le rayon de courbure de l’hypercirconférence en un point (uw, v) est donné 
par la formule: 


m u \3/2 
1 Ia 
0 ] ( ape 1? a) 
u m2 uU 
4 
ch 4 P sh an 


Aux extrémités du diamétre vertical (w = 0) le rayon de courbure de l’hyper- : 
circonférence est égal a —/, donc indépendant du module, tandis qu’aux 
extrémités du diamétre horizontal le rayon de courbure est égal a 


[2 
a eee ou j 22 


en appelant y, l’ordonnée canonique de la corde horizontale qui découpe une 
caténe de longueur 2/ dans la chainette de module m. Ce rayon de courbure 
dépend donc du module. Cette différence remarquable entre ces deux rayons 
de courbure extrémes s’explique géométriquement. Un diamétre vertical cor-* 
respond en effet a un arc de chainette de longueur /, situé a l’infini sur la 
chainette et dont la corde a la méme longueur que celle de l’arc. Une corde 
infiniment voisine procéde de la précédente par deux translations des extré- 
mités, translations qui sont égales et sensiblement confondues avec l’arc, donc 
avec sa corde. : 

D’autre part, dans la position a l’infini, la corde est normale a la bissectrice 


des deux tangentes a ses extrémités, donc les deux rayons vecteurs OP normaux 
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a la tangente sont égaux. Il n’en est pas de mérie avec les cordes horizontales, 
ot les cordes infiniment voisines sont bien normales a la tangente mais ne sont 
plus égales parce que les translations de leurs extrémités bien qu’également 
inclinées sur l’horizontale ne sont plus confondues avec celle-ci et de plus sont 
inclinées en sens contraires. Cette inclinaison influe donc sur la longueur du 
rayon de courbure, ce qu’exprime l’apparition du terme y,/m dans la formule 
du rayon de courbure a une extrémité du diamétre horizontal. 


Figure 5 


Par un raisonnement géométrique on établit une autre expression du rayon 
de courbure d’une hypercirconférence ou apparaissent les rayons de courbure 
a et b de la chainette aux deux extrémités A et B de la caténe qui forme 
Vhyperrayon, et l’angle 9 des deux tangentes a la chainette a ces extrémités. 


Cette expression est 
Oo mons 
ee Toe ee 


g=4 


Cette propriété fondamentale de la tangente en un point d’une hypercircon- 
férence incite A étudier comment se répartissent sur une chainette toutes les 


extrémités A et B des arcs AB, dont les tangentes aux extrémités forment des 
angles dont la bissectrice a une direction invariable 1 (figure ay. 

Soit alors J le point de la chainette dont la tangente a cette inclinaison #, et 
soit J le point ot la tangente a l’inclinaison i+ hi j (nous dirons que ces 
points J et J sont en correspondance orthogonale). On démontre alors que les 
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extrémités A et B d’un arc dont la bissectrice tangentielle (bissectrice de 
angle formé par les deux tangentes a l’extrémité de l’arc) a la direction 7 
sont liées par la relation: ads 

fB IB 

fA [A 

Sur la chainette les deux points A et B sont donc en correspondance an- 
harmonique avec les points I et J (qui eux sont en correspondance ortho- 
gonale) définissant la direction de la bissectrice. 

On retrouve ici une analogie avec la correspondance anharmonique qui 
existe entre les points d’intersection sur la base commune a tous les triangles 
de sommet donné et constitués par deux autres cdtés issus de ce sommet et 
formant des angles dont les bissectrices ont une direction fixe. 

Par un point extérieur 4 une hypercirconférence on peut toujours mener 
deux hyperdroites, équimodulées, qui lui soient tangentes. Si cette hyper- 
circonférence a un hyperrayon dont la longueur est supérieure 4 son module, il 
existe alors aux extrémités de son diamétre vertical une zone intérieure, par 
les points de laquelle on peut tracer une hyperdroite tangente a l’hypercircon- 
férence. 


Normale et hypernormale 


La normale a une courbe en un point de celle-ci est le lieu des centres de 
toutes les circonférences qui sont tangentes a cette courbe en ce point; l’hyper- 
normale a une hypercourbe en un point de celle-ci sera le lieu des centres de 
toutes les hypercirconférences d’un module donné mais d’hyperrayon quel- 
conque et qui sont tangentes a l’hypercourbe au point considéré. 

Si on remarque que dans le cas de la normale a une droite, les circonférences 
peuvent étre assimilées a des hypercirconférences de module infini, donc ayant 
méme valeur que celui de la courbe (la droite) a laquelle elles sont tangentes, 
nous pourrons préciser que l’hypernormale a une hyperdroite est le lieu des 
centres des hypercirconférences d’hyperrayon quelconque mais de module égal 
a celui de l’hyperdroite (chainette), et tangentes a celle-ci en un point donné. 

D’ailleurs la remarque faite im fine du paragraphe consacré aux hyper- 
circonférences permet déja de conclure que si on fait varier le module d’une 
hypernormale (c’est-a-dire le module des hypercirconférences tangentes) on — 
passe de l’une a l’autre par une homothétie dans le rapport des modules ayant 
pour centre le point de contact. 

Cette déformation homothétique des hypernormales, quand leur module 
varie, fait passer progressivement celles-ci de la droite verticale (module zéro 
des hypercirconférences) 4 la droite orthogonale a la tangente a la chainette 
au point considéré (module infini des hypercirconférences), les formes inter- 
mediaires n’étant nullement des droites. 
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Sous le bénéfice de cette remarque nous étudierons l’hypernormale équi- 
modulée, c’est-a-dire celle qui est le lieu des centres des hypercirconférences de 
module égal a celui de ’hyperdroite. 

Une hypernormale équimodulée 4 une hyperdroite ne peut couper celle-ci: 
il faudrait pour cela qu’il existat sur une chainette une caténe dont la bissectrice 
tangentielle (bissectrice de l’angle formé par les deux tangentes aux extré- 
mités) soit paralléle 4 l’une des tangentes d’extrémité, ce qui est impossible. 
Mais une hypernormale équimodulée 4 une hyperdroite devient tangente a 
celle-ci au point a l’infini sur la branche opposée a celle ot se trouve le pied 
de l’hypernormale. Si l’hypernormale est tracée sous un module supérieur a 
celui de l’hyperdroite, alors elle coupe celle-ci; inversement, tracée sous un 
module inférieur elle ne la coupe pas, mais lui devient paralléle — et nontan- 
gente — au point a l’infini (figure 6). 

Soit x l’abscisse canonique d’un point d’une hyperdroite (chainette) de 
module m. 

Rapportée a des axes orthogonaux, paralléles aux axes canoniques de la 
chainette et ayant leur origine en ce point, l’hypernormale, équimodulée, a la 
chainette en ce point a pour équation: 


Rapportées aux axes canoniques, les coordonnées 
x 
y Spm ch, X=uUu+% 


de cette hypernormale satisfont donc a l’équation: 


AS 
ai oh cha Wagienbulinl 5 Teebel 
m m 


RIES wien dy 
m 


Toutes les hypernormales aux divers points d’une hyperdroite, équimodulées 
avec celle-ci, ont une enveloppe dont l’équation est: 


Y=mch> +m [1 + (nye 


Entre l’abscisse canonique X du point caractéristique dune hypernormale 
4 une chainette et l’abscisse x du pied de cette hypernormale on a la relation: 


ae tog es 
CL rata sh ae 
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= 


ated 


Figure 6 
Hypernormale de module m, élevée 4 une hyperdroite de module m, en un point P de celle-ci. 


a5 
Equation de l’hyperdroite y=m Cee R 


u 
—m sh — 
m % £4 
équation de l’hypernormale v= 3; U = margsh [n= on | di 


& ™m 
sh— 
mm 


Le point caractéristique d’une hypernormale marque le centre d’une 
hypercirconférence dont le point de tangence avec la chainette (le pied de 
Vhypernormale abaissée sur l’hyperdroite) est de l’autre cété de l’axe vertical, 
Vhypercirconférence et la chainette ayant en ce point méme rayon de courbure. 
Cette égalité des rayons de courbure se démontre géométriquement, et peut 
aussi étre montrée par le calcul. : 

La longueur / de l’hyperrayon de l’hypercirconférence tangentielle ayant 
pour centre un point (X, Y) de la courbe enveloppe des hypernormales est donnée 
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par la formule: 
22 -(Y —™ ch \* + 2 ir (ny fv Hi ch)? = m?| , 


La courbe enveloppe des hypernormales A une hyperdroite partage le plan 
en deux zones. De tout point de la zone supérieure — que nous appellerons zone 
tripode — on peut tracer trois hypernormales a l’hyperdroite, savoir une dont 
le pied P, est du méme cété de l’axe que le point, et les deux autres ayant 
leurs pieds P, et P; de l’autre cété de l’axe (figure 7). 

Ce point tripode est donc le centre de trois hypercirconférences tangentes 
a V’hyperdroite. Celle dont le point de tangence est du méme coté de l’axe a 
Vhyperrayon le plus court et le point de tangence le plus élevé au-dessus de la 
base. Celle dont le point de tangence est intermédiaire entre les deux autres est 

C 


Py 


| 
I 
gp yeep ot See pee ae 
Figure 7 
Point tripode A. Phe wet 
i i téristiques C1, Cy, Cg. Ses yper- 
Ses 3 hypernormales, de pieds P,, Pz, Pg et de points carac ues Cy 
P acallates de points ae Per oecemaent R,, Ry, Rj (onn’a figuré que la moitié de chaque hyperparalléle). 


\ 
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celle dont I’hyperrayon est le plus long, et dont le point de tangence est aussi 
le plus bas; son rayon de courbure en ce point est supérieur a celui de la 
chainette, ce qui n’est pas le cas pour les deux autres hypercirconférences con- 
centriques tangentes a la chainette. En résumé des trois hypercirconférences 
tangentielles de centre N, l’une, celle 4 point de tangence P,, est entiérement 
dans la zone concave de l’hyperdroite (chainette), celle 4 point de tangence F, 
est a l’extérieur de l’hyperdroite dans sa partie inférieure et ne pénétre dans la 
zone concave qu’en des points situés au-dessus de P, et de P;; enfin Uhyper- 
circonférence a point de tangence P, est a l’intérieur de la zone concave aux 
abords du point P;, mais coupe l’hyperdroite en un point compris entre P, 
et P, pour sortir de la zone concave et y rentrer en un point de l’hyperdroite sis 
au-dessus de P,. 

D’un point de la zone intérieure — que nous appellerons zone monopode — 
on ne peut mener qu’une seule hypernormale a l’hyperdroite. 

D’un point de la courbe enveloppe on peut évidemment mener deux hyper- 
normales, l’une courte ayant son pied du méme cété de l’axe que le point, 
tandis que l’autre longue a son pied de l’autre cété: le rayon de courbure de © 
Vhypercirconférence tangente a la chainette correspondant a l’hypernormale 
longue est, au point de tangence, égal a celui de la chainette. . 

Ces particularités des hypernormales sont analogues a celles que présentent 
les normales a une chainette, avec cette difference que dans ce dernier cas la 
courbe enveloppe coupe la chainette, étendant ainsi la zone tripode a l’extérieur ~ 
de la zone concave de la chainette. 


Paralléles et hyperparalléles 


Utilisant la méme analogie que celle qui nous a fait passer de la définition 
de la normale a celle de l’hypernormale, nous définissons l’hyperparalléle 4 une 
hyperdroite tracée A une hyperdistance / équimodulée, de celle-ci, comme le | 
lieu des centres de toutes les hypercirconférences, équimodulées également, 
tangentes a cette hyperdroite et dont les hyperrayons ont une longueur com- 
mune J, 

C’est en somme la ligne que décrit le centre d’une hypercirconférence quand 
on fait glisser celle-ci tangentiellement (en la maintenant verticale évidemment) 
le long d’une chainette de méme module (figure 8). ; 

Il est aisé de démontrer qu’en tout point d’une hyperparalléle la tangente | 
est paralléle a celle a la chainette au point ot l’hypercirconférence qui admet le © 
premier point pour centre est tangente a cette chainette. 

En effet si on considére deux centres voisins Cy, C, de deux hypercirconfé- | 
rences égales et tangentes lune en P,, l’autre en P, a l’hyperdroite, on peut | 


dans l’hypercirconférence de centre C, tracer le vecteur Cy fi, €quipollent au 
vecteur C0). 
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D’autre part on a l’égalité géométrique: 


Co Cy = Co t+ Py Pyt+ PC, =P, P+ C, Py +h, C= P,P +h, 


et a la limite quand C, est trés voisin de Cy, les deux vecteurs P, P, et $1 Py 
étant confondus en direction, avec la tangente commune en P, a la chainette 


Ia Br Saeco 4 


Figure 8 


Hyperparalléle, 4 une hyperdistance /, 4 une hyperdroite de module m (cas de //m = 8). 


emplese 


2 


paralléle a cette tangente. 

De cette égalité géométrique, combinée avec le parallélisme des tangentes, on 
. déduit que le rayon de courbure 4a l’hyperparalléle en un point Cy est égal a la 
différence des rayons de courbure en Pj, ala chainette et a l’hypercirconférence. 

Tant que l’hyperdistance /, qui définit le parallélisme, est moindre que le 
module m de l’hyperdroite, les hyperparalléles correspondantes présentent 
lV’aspect d’une ligne continue, ne coupant qu’une seule fois l'axe de l’hyper- 
droite, avec une tangente horizontale en ce point. 

Pour les hyperdistances supérieures au module m, la forme des hyperparal- 
léles est celle d’une ligne qui coupe trois fois l’axe de l’hyperdroite. La premiére 
- intersection, celle-la 4 tangente horizontale, marque un point a une distance / 
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au-dessus du sommet: les deux autres intersections, croisées, se font 4 un point 
supérieur. 

Ces hyperparalléles se terminent vers le haut par des branches qui tendent 
a devenir paralléles a l’hyperdroite, et a une distance de celle-ci égale a 


l 
2m argsh = 


(longueur du rayon vecteur horizontal). 

De part et d’autre de l’axe et a une cote inférieure a celle du franchissement 
horizontal de l’axe, l’hyperparalléle présente un point de rebroussement, sis 
évidemment sur la courbe enveloppe des hypernormales: |’hypercirconférence, 
dont le point de rebroussement est le centre, a en son point de tangence avec 
la chainette méme rayon de courbure que celle-ci en ce méme point. Le rayon 
de courbure y est bien nul comme nous |’avons démontré. 

Et on aboutit a cette conclusion que par un point situé dans la zone tripode 
des hypernormales on peut mener trois hyperparalléles 4 une hyperdroite, et 
deux seulement si le point est sur la ligne enveloppe des hypernormales; et une, 
si le point est dans la zone monopode. 

Les coordonnées Y et X d’une hyperparalléle 4 la distance / sont données, 
en fonction des abscisses canoniques x des points de la chainette, par les for- 
mules: 

A= x + margsh(— sh), Y = mch— +0 
m m m 


ou v est exprimé par 


v2 = [24-2 m? eh at |/m* cht a [2 sh2—. 
m WM m 


Ces formules permettent de retrouver, au terme de laborieux calculs, la 
propriété fondamentale du parallélisme des tangentes, aisément démontrée par 
la géométrie. 

La figure ci-contre montre la forme générale d’une hyperparalléle A une 
hyperdroite, tracée a une hyperdistance de celle-ci supérieure 4 son module m 
(figures 8 et 9). 


Pour résumer cette étude de l’hyperparallélisme nous remarquerons que - 


pour des étres raisonnables dont les perceptions visuelles ne dépassent pas le 
module de leurs régles (instruments) de mesure il est conforme a leur bon sens 
que par un point on ne puisse mener qu’une seule paralléle 4 une droite; ce 
n’est que quand la portée de leurs perceptions dépassera ce module qu’un 
nouveau bon sens leur fera trouver tout naturel que par un point on puisse 
mener trois paralléles — et a des distances inégales — a la ligne qui leur tient lieu 
de ligne droite, c’est-a-dire de tracé de mesure. 


——— 
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Hi(2, > m)| 


Hii, < m) 


tay 
| 
| 
| 


Sie cel 


Figure 9 


(Hil, < m): Hyperparalléle 4 une hyperdroite tracée 4 une hyperdistance Ig<m. 
(Son hyperparallélisme différe peu en apparence du parallélisme euclidien.) 


Cas particulier de la courbe enveloppe ®, des chainettes isosthénes 
d’un faisceau 


Le role fondamental tenu par cette courbe enveloppe dans |’étude des pro- 
blémes mécaniques posés par le fonctionnement des téléphériques et des 
télésiéges, nécessite un examen spécial quant aux hypernormales et hyper- 
paralléles 4 cette courbe. 

On sait que chaque point de cette courbe marque le contact avec une 
chainette isosthéne du faisceau, dont le module m est d’autant plus grand que 
la cote du point au-dessus de la base du faisceau est plus élevée, en méme 
temps qu’est plus voisine de la verticale la tangente (i, désignant l’angle que 
celle-ci fait avec l/horizontale (figure 10). 

Ainsi a la différence d’une hyperdroite, il convient d’assigner 4 chaque 
- point @; de la courbe ® un module particulier, celui (m,) de la caténe isosthéne, 
dont ce point est le point caractéristique. 

A chaque point de la courbe ® correspondra donc une hypernormale 
définie par son pied y,, son module m,, et l’inclinaison 7, de l’élément tangentiel 
auquel doivent étre tangentes les hypercirconférences de module m,, dont les 
centres traceront cette hypernormale. 

Ces hypernormales (y,), dismodulées (c’est-a-dire dont le module varie 
avec leur pied g;) ont une courbe enveloppe, qui est le lieu de leurs points 
caractéristiques; le point caractéristique C, de l’hypernormale (y,) est le 
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centre de l’hypercirconférence de module m,, tangente en gy, a la courbe @ 
et y ayant un rayon de courbure @,, égal a celui de la courbe @ en ce point. 

La connaissance des rayons de courbure de la courbe est donc importante; 
on les calcule aisément en deux points particuliers, savoir au sommet ou il est 
égal A la hauteur / du faisceau et au niveau h, ot il est égal a 1,810 h. 


ra) 


Figure 10 
® courbe enveloppe des chainettes isosthénes passant par A (A a =h). 


Hyperellipses 


Soient deux points /, et & et L une longueur. Une caténe de longueur L, 
tracée de /, a & aura un module m; bien déterminé qui sera fonction de la 


longueur du vecteur rectiligne /, & et de l’inclinaison de celui-ci au-dessus de 
Vhorizontale. 


_—-~ 


Tous les points de cette caténe /, fy et aussi ceux de sa symétrique inversée 
(symétrique par rapport au milieu M, de F, Fj) sont des points dont la somme 
des hyperdistances /, A + , A, mesurées sous module m, est égale a L. Ces 


deux caténes, [4 /, et son inverse, figurent donc l’hyperellipse, d’hyperlongueur 
L, en module m. 
On démontre que si l’on fait croitre le module m, le lieu des points dont la - 


— 


eet 
somme des hyperdistances F, A + F, A mesurées sous ce nouveau module est 
égale 4 L, se compose de deux lignes continues, symétriques par rapport au 
milieu M, et enchassées les unes dans les autres, celles de plus grand module 
entourant celles de module moindre, la forme limite extréme étant 1’ellipse 
classique correspondant a des hyperdistances mesurées sous un module infini. 


— 


tai aes} 
Si au contraire on mesure les hyperdistances F, A + K, A sous un module 
moindre que m, on démontre que ce module ne peut étre moindre que m/2, 
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auquel cas l’hyperellipse se réduit au point M, milieu de i, fy, tandis que pour 
les modules intermédiaires entre m et m/2 les hyperellipses se présentent 


comme des courbes progressivement enchdssées les unes dans les autres (fi- 
gure 11). 


- <m< m,) (m,) (m™ 7 oa 


(m=) 


‘ Figure 11 
Hyperellipses de foyers Fj, F, et longueur L = FA + FA . 


a 


FA +F,A = L|(m,); FA + FA =L](m =00); MF, + MF, = L|(m)2). 


Summary 


Differential calculus is proving inefficient to solve problems of mechanics 
dealing with ropes, unless one satisfies oneself with rough approximations. 

Ordinary or Euclidean geometry is also proving useless: The author shows 
that in order to solve such problems as described before, one has to use a special 
branch of geometry. Distances are no longer measured along rigid straight rules, 
but along the catenary curves in question. 

The fundamental characteristics of this hypergeometry, with some of their 
more striking results, are described in the following study. 


(Regu: le 11 aotit 1956.) 
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Stresses Due to a Nucleus of Thermo-Elastic Strain 
(i) in an Infinite Elastic Solid With Spherical Cavity and 
(ii) in a Solid Elastic Sphere 


By Brauma Dev SuHarma, Pilani, Rajasthan, India?) 


Abstract 


In this paper solutions in series form for the stresses due to a nucleus of 
thermo-elastic strain in an infinite elastic solid in the presence of a spherical 
cavity and also in an elastic solid sphere have been found. 


Nomenclature 


The following nomenclature has been used in this paper: 
x, ¥,2 Cartesian coordinates; _ 
v,0,q_ spherical polar coordinates; 
Up, Uy, U, Components of displacement in Cartesian coordinates; 
U,, Ug, U4, components of displacement in spherical coordinates; 
6; 06, Sy, Tro, Trg Top Components of stress in spherical coordinates ; 


E coefficient of elasticity in stress; 
G coefficient of elasticity in shear; 
a coefficient of linear expansion ; 
vy  Poisson’s ratio. 


Introduction 


The thermal stress problem of an infinite elastic solid at zero temperature 
except for a heated region has been solved by GooptER [1]?). In this paper the 
infinite solid is considered to have a spherical cavity while the heated element 
is at some finite distance from it in the solid. Numerical results for the stress ~ 
and displacements on the surface of the cavity are given. The results have been 
compared with the corresponding numerical values appropriate to the two- 
dimensional analog of the present problem. A second problem of a solid sphere 
at zero temperature having a heated nucleus inside has also been solved. The 
stresses and displacements on its surface have been calculated. 


1) Birla College. 
*) Numbers in brackets refer to References, page 150. 
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Solution 


Consider an infinite elastic solid at zero temperature, while an element of 
volume dQ at (c, 0, 0) is heated to a temperature T. The displacement in this 
case 1s given by the gradient of y where 


ke Si 1 
ae ge (r2 — 2recos6 + c®)® ? (1) 
in which 
A=-2G oe and m= ea 
4% —y 
Hence 
2Gu, = MA ri: 005:0) nis | 
(x? — 2recos8.+ c)32 ” 
20a é Acsin@ (2) 
(v2 — 2rccosé + c?)32 ? 
u, =O, 
and 
ah eer 272+ 3c* cos? @— 4rccos6 — c? 
oi (yr? — 2rccosé + c?)>2 4 
se a 3 c? sin? 6 ‘ 1 
6 * |? = 27 c cos6 + c?)F? EXE eat 
s 1 3 
a A (v2 — 27 ccos6 + c?) ” 3) 
A 3c sin 6 (y — c cos 6) 
Trg = 


ae 7 — "92 £:C-c0si0 = c7)P!2 © 
Tre = Tp, = 0. 


We denote this displacement and stress field in the infinite solid by [y]. 

The stress in the infinite solid being known, in both the problems under 
consideration we have to free the spherical surface of traction. In the first 
case the nucleus is outside the surface while in the second case it is inside it. 
These problems can be attacked by any one of the various stress-function 
approaches to the axisymmetric problems of elasticity. We shall utilize the 
approach originated by BousstnesQ according to which the general solution 
of displacement equation of equilibrium, in case of torsion free rotational 
symmetry, and the absence of body forces, is representable as the sum of two 
displacement fields. Ae oe a ees (4) 
2G (Ug, Uy, Me) = grad (zy) — 4 (1 —») (0,0, 9), (5) 


where V2 9 =V? y=0. 
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These basic displacement and associated fields of stress can be transformed 
into general orthogonal axisymmetric coordinates. We here obtain the parti- 
cular solutions by introducing the interior and exterior spherical harmonics of 
integral order as generating stress functions m and y [2]. 

We consider two pairs of aggregates of harmonic stress functions 


g(r, 6) =7-"—" P. (cos6) ,, “Yj by—=7 © Poicos0); (6) 
and 
yp, (r, 0) =” P, (cos6), y,(r,0) =” P, (cos8) , (7) 


in which P,, denotes the Legendre polynomial of order . The stress functions 
defined in (6) represent exterior spherical harmonics while those in (7) represent 
the interior ones. The particular solutions generated by @y, Wn, Yn and yy 
will be distinguished by [A,], [Z,],[C,] and [F,] respectively. Considerable 
simplification arises from a replacement of solutions [£,] and [F,] with the 
linear combinations 


[B,] = (2m + 1) [£,]—(#+4—49) [4,1], (8) 


[D= 20> 1) lia 3p iC pee (9) 


We now give the displacement and stress fields appropriate to the solutions 
[A,.], [By], [C,] and [D,]. With the auxiliary notation 


) = Cos 0, b= sin 6 (10) 


we obtain the solutions for [A,,] 


1 ie , 
2Gu,= Ph, 2.G y= — Loe, Uy =0, (Lh 
_ (m+ 1) (n+ 2) 
a ynt+3 Pas | 
1 ’ 
Oy = sare (Pa ~ Fd) (9 4 2) PI 
P’ 
eet (12) 
Vi 2. ae, 
to = mee PE,» 


Upp = tog = 0% 
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for [B,,| 
a cle ty (n+ 4 $y 
264, yrtl ap eae 
,n—3+ 4 , 
2G Uy = P= ynth * i% tl? (13) 
u,=O0, 
0; = (n +1) [(m +1) (n +4) — 20] 2242 
5 yn+2 
1 
Bor pare tA) (et od 2) Pie (3A Py), 
pee ee: fang 1)(2n+1)P = (14) 
Oy Sena, v) (n+1)(2n+1)P.,4+(n—-34+4»)P), 
Tro = (n?-+ 2n—1+ 29) p—H2, 
Trp =, = 0, | 
for [C,] 
264,—=nr" P) 26u,=—-pr'P,, 4u=9, (15) 
6,=n(#—1)7-"*P,, | 
6, = 9" (Po — (+0 +1) PY, 
Oy _ kee [—Pies + (2 n ie 1) 2 ie , (16) 
Tro =—(n— 1) pr Pr, 
Tg = Tre = 9, 
for [D,] 
264, = nw(n—3 +49) 4% Piog; 
2G y= —(n+4—40)"PL,, (17) 
u,=O0, 
o, = [n? (n — 3) —2yn]r"P,_,, 
j= [(n+4—4) Pi —n(n? + 3n+3—2) Fy), | 
p= [n+ 4-49) P+n(2n+1) (1-2) Pil, ¢ (18) 


P 


9 = (2—? — 24) 14 Pl, 


Tyo = Tp = 9, 


_.. ZAMP VIII/10 
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Infinite Elastic Solid With Spherical Cavity 


We take the centre of the cavity as the origin and the radius to be unity. 
The heated particle is at (c, 0, 0) such that c > 1. In order to free the cavity 
from traction 


Cn", et, 0 -¥ Sa, 0 SU= a Oey =a (19) 


Also, all the stress components must be zero at infinity. 

In the region under consideration, 7 = 1, the stress systems [A,,] and [B,] 
are uniform, and hence we can superpose them on the solution [y] for the 
infinite solid. We assume solution [S] to the problem characterized by the 
boundary condition (19) in the form 


[S]=xt D7 { 4l4,] + 2,[B, 1] }- (20) 


Now, for the solution [y] on 7 = 1 


2 + 3 c2 cos?0 — 4¢ cos — c? 
o,=/,(0) =A (1 = 2¢ cos® + c?)P?2 ‘ 


; (21) 
3c sin6@ (1 — c cos6@) 
(1 — 2c cos6 + c?)b/2 © 


T6@ = fo(9) =A 


f,(0) and 7,(0) can be expressed in Legendre series and when c > 1 we get: 


= iP. i) 
(0) =A 3 mln — 1) =e 
n=1 
(22) 
— @& pe fa) 
f(0) =A pS’ (w—1) Sa eos) 
n=1 
A little reduction shows that condition (19) is satisfied if 
a, (n+ 1) (n+ 2)+6, 0 (nw? +3n—2r)+A AE) =); 
an (+ 2) +b, (2-2 420) 4M=) _o, 
or 
Ripe. A n(n —1) (2n—1) 
ey Z{1ltn+ n7—9(2n+1)}crt1? 
Bites A (n—1) (2n+1) (23) 
2{l+ n+ n?—»(2n+1)}crt1 
REOPEN 


The value of these coefficients together with equations (20), (2), (3), and (11) 
to (14) constitute the complete solution of the problem. 
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Now, we find the numerical results when » = 1 /4 and c = 5. The displace- 
ment components on the cavity are shown in Figure 1. Figure 2 shows the 
stress concentration on the cavity. There are two rings on the cavity where og 
is zero. It is interesting to note that the maximum value of 0, — 09 Occurs near 
about the same value of 6, 70° nearly, where wu, vanishes. The stresses around a 


Figure 1 
Deformation of the spherical cavity. 


0-06 


0-04 


-0-02 


Figure 2 
Stresses on the surface of the spherical cavity. 
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circular hole in an infinite plate are compared with those in the present three- 
dimensional problem in Figure 3. As usual, the stress concentration in space 
problem is found to be of less severe intensity and more localized in character 
than it is in its two-dimensional counterpart [3]. 


0-06 


INFINITE SOLID 


INFINITE PLATE 


Figure 3 


Comparison of the stress on the spherical cavity in an infinite solid with that on the circular hole 
in an infinite plate. 


Solid Elastic Sphere 


As in the first problem, we take the centre of the sphere as the origin and 
the radius to be unity. The heated nucleus is at (c, 0, 0) such that ¢ < 1. In 
order to free the surface of traction, again 


Oy eG eg Oe Oa OS rie rer Zien (24) 


Inside and on the surface of the sphere, the stress systems [C,,] and [D,] are - 
regular and we can superpose them on the solution [y] for the infinite solid. 
We assume solution [S’] to the problem characterized by the boundary con- 
ditions (24) in the form 


For the solution [y], 6, and og on 7 = 1 are given by (21). In the present problem 
c < 1, and hence the representation of /,(0) and /,(0) in Legendre series in this 
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case will be 


(Olax = f(0) = A 3 (n+ 1) (WF 2) 0 P, (0080), | 


n=0 
= (26) 
[olen = fo(0) =A p D” (m + 2) c* Py (cosé) . | 
n=0. 
It can be easily seen that (24) is satisfied if 
pee Fe a 8) n 
aS CE 6 es “SEEN es ; 
eee (m + 2) (2 +1) hs 27 
dy, Mee ati cseeri ( ) 
(= Of, 2... ops 


With these values of the coefficients c,, and d,, equations (25), (2), (3) and (15) 
to (18) constitute the complete solution of the problem, 

Here also the numerical results for y = 1/4 and c = 1/5 have been calculated. 
c = 1/5 is reciprocal point to that in the infinite solid problem with respect 
to the sphere of unit radius. The displacement and stress components on the 
surface are represented in Figures 4 and 5 respectively. It is interesting to 


Figure 4 


Deformation of the surface of the solid sphere. 


note that in this case #» is much smaller than wu, on the surtace, but c up is 
of comparable size. The value of o — o, is much smaller in comparision to 
Oo and oy. 
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Figure 5 


Stresses in the surface of the solid sphere. 
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Zusammenfassung 


Die thermischen Spannungen in einem festen K6rper unendlicher Ausdeh-.- 
nung, welcher einen spharischen Hohlraum enthalt, sind bei einer Temperatur 
von 0°C in Gegenwart eines erhitzten Elementes, das sich in endlichem Abstand 
vom Hohlraum befindet, hergeleitet worden, wobei zahlenmassige Angaben fiir 
die Spannungen und Verschiebungen an der Oberflache des Hohlraums gemacht 
werden konnen. Die Ergebnisse sind mit den entsprechenden, fiir den zweidimen- 
sionalen Fall giiltigen Zahlenwerten verglichen worden. Ferner war es méglich, 
auch fiir das Problem einer festen Kugel: von der Temperatur 0°C und einem 
erhitzten Kern in ihrem Innern eine Lésung zu finden. 


(Received: Juni 6, 1956.) 
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On a Paper of Gaunt Concerned With the Start of Numerical Solutions 
of Differential Equations 


By RENE DEVOGELAERE, Notre Dame, Indiana, USA!)2) 


Introduction 


In 1927 Gaunt [3]*) published a very important result which seems not to 
have attracted any attention. He showed that if the difference method 


be 0§; = h f(t, é),, = Xo) &, 
is used to solve the differential equation 
= f(t, x), (0) = 4%, 


the value of £,, which is necessary to start the method, should be computed in a 
specific manner not coinciding with the exact value of ¥(#). The error will then 
behave uniformly in the sense given below and will not have the oscillatory 
character which is often introduced if &, = x(h) is determined by the Taylor 
expansion of ¥ at the origin. Moreover, it then becomes possible to use RICHARD- 
sON’s metamethod which permits to obtain a better approximation from two or 
more computations with different intervals h. 

I would like to state his result in a precise form and give a more elegant proof. 

I will also state a similar theorem for second order differential equations. 
Gaunt’s result can be generalized to many other difference methods to solve 
differential equations ordinary and partial. 


1. A Theorem of Gaunt 


Theorem I. /} f(t, ~) has continuous second partial derivatives, in a domain D, 

the solution x of the differential equation 
a = f(t, x) (1.1) 
with initial condition (ty, X)) © D is approximated in the domain D by the solution 

of the difference equation 

wOE, =hflt, &), teat +th, S=%, (1.2) 
h2 
a 


=a that aH +o%m+¥e,, lim lel=0, (1.3) 


in such a way that 
E(t) = a(t) +h? (nt) + oh A] (1.4) 


where y(t) is a continuous function and C(t,, h) tends to zero with h. 
1) Department of Mathematics, University of Notre Dame. 


2) Under contract with the Office of Naval Research. 
3) Numbers in brackets refer to References, page 155. 
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More specifically y(t) is the solution of the differential equation 


; 0 i awn 4 
n= ae hier wv, nt) =O; (1.5) 


where x must be replaced by the above solution of (1.1). 

I will give the proof without using the known relation between the operator 
uo and the operator D of differentiation with respect to ¢ and will indicate in the 
next theorem, which asks for stronger conditions, how to use this relation to 
generalize (1.3) and (1.5). 

Under the hypothesis of the theorem, it is clear that 7 is continuous, hence 
that the solution of (1.5) exists, is unique and has a continuous derivative. 
Hence we have the Taylor expansions 


: nae eae 
ai + h)= «(i) + halt) + = #() + [4 @) + el, | 


y(t + h) = n(t) + A[y(t) + e], _ (1.6) 
ft, w+ Ax) =f, + Ax( oF i e) ; | 


where « tends to zero with h. 
(1.6), for instance, permits us to write, with +; = +(t,) , 


p O%, = Myr — r= he + (HF 8) - (1.7) 


Hence, if we define C(é, h) by (1.4), replace in (1.2) & by (1.4) and use (1.7) 
as well as similar relations, we have 


F hes h 


af (1.8) 
Because of (1.1) and (1.5), ¢; satisfies the relation 
h 0 
Maes seth |(S) + 6] t+ hem, (1.9) 


with ¢, = Oand |¢,|< h (e, + 2 e); indeed, (1.3) and (1.5) at ¢ =#, give: 
h2 h2 h3 


h? by = 4) + h i, 4 : Ky + h8 ey | 4 hat 3 ate 6 (4% + €) 


—h® [ng +h (ing + &))- 


Therefore, it is possible to prove that ¢; tends to zero with h, and uniformly so 
if the solution of (1.1) remains in the domain D. 


2. Generalization of Gaunt’s Theorem 


The theorem given above does not cover exactly the result of Gaunt but is — 


contained in it, the result of Gaunt is closely related with the theorem given 
below, when ” = 3. 
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Theorem II. If f(t, x) has continuous n-th partial derivatives in a domain 
D,n> 2, the conclusions of Theorem I hold with (1.3) and (1.4), replaced by 


h2 
= %+hx + Ar Rik a «eee eal lim e,= 0, (2.1) 
1 Iv eee Of 
=——— |x, + 2 (Ge)] >. 2.2 
Ya 24 : %o OX 2x4) 
and 
E(t) = x(t) + B® ny(t) + AA g(t) +e + 12? ny (t) + AN LU, A), (2.3) 


where 2p is the largest even integer smaller or equal to n and €(t,, h) tends to zero 
with h. 


For the proof of Theorem II, it is easier to use the difference operators to 
prove that the odd powers of # disappear in (2.3) and to generalize (1.3) and 
(1.5). I will do so but will not consider the remainder terms which are to be 
treated exactly as in Theorem I. I will recall the formula 


h3 D3 pent) pent 


po=sinhh D=hD > A +t + (2n + 1)! tee, (2.4) 


If we replace in (1.2) € by its value (2.3) and use (2.4) we obtain, after 
division by h, (7) = *) 


i p2t D+ of [ ; 1 oz [ Tae 33 
Dr eg ary 8 t ae | OME te) +3 x Ge On) 


ad t= ¢=1 
(2.51) 
This gives the equation (1.5) for 7,, and 
a, RE 4) 
Ese oe sigh Ac ae iat mi Br? No(%) = 0. (2,52) 


We then prove by recurrence that (2.3) holds. We only need to remark that 
(2.3) is correct when ¢ = 0, which is obvious and when ¢ = —h which will deter- 
mine 7_,, hence 7,: 


62 
f= fo + w 08) + —- So 


or 
Lae t+4) ses 
& = Hy +h X%o+ 5 » ae Tee eo 


t=17= 


because (2. 3) is satisfied for = 0 and t= h and must be satisfied for ¢= —h and 
because of 


hei af 
Ce es Aaa 
This gives (2.1) and (2.2). 
3. Comments 


The importance of the contribution of GAUNT will become clear from the 
following comments. 
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(i) The above method is not restricted to that method of difference considered 
by Gaunt; it shows that contrary to the general belief (and our own belief 
before we read this paper), one cannot start a difference method by any pro- 
cedure: an other difference method, Taylor expansion, Runge Kutta method. 
Indeed, each difference method will ask for the determination of initial values 
according to well defined formulae which do not necessarily coincide with the 
exact value of the solution at the corresponding point; we see, for instance that 
(1.3) differs from the Taylor expansion of *(/) at the third order, hence the use 
of Taylor expansion would give, in general, an error of the form h? C(¢, h), where 
¢ is bounded. On the contrary by using the correct procedure the error is made 
uniform in the sense of formula (1.4) or (2.3). An example is given by GAuNT. 

(ii) Theorem II does prove that the difference method described in this 
theorem is of class RC n,{2,4,..., 2p} (see[1]), hence that RICHARDSON’s 
metamethod can be applied [1], [5]. An example is given by Gaunt [3]. 

(iii) The choice of the difference method made by Gaunt was a lucky one, it 
was guided by an unconscious realization of the importance of the symmetry 
of the method. (Central difference methods are most of the time symmetrical 
by their very nature and it is this symmetry which leads to error formulae like 
(2.3) where all odd powers of f are missing.) The advantage of the form of (2. 3) 
is especially great when RIcHARDSON’s metamethod is used. This point is discussed 
in detail in [1]. 

(iv) The above form of solving differential equations with difference methods 
plus RicHaRDSON’s metamethod constitutes an excellent way of checking com- 
putations, especially precious when computations are not part of a large 
program (see [2]). 

(v) Gaunt’s method constitutes an excellent way of studying the errors in 
the large of a method. 


4. Computation of the Starting Value 


Clearly the use of formulae (1.3) and (2.1) presents the same computational 
difficulties than the use of Taylor expansion except perhaps that in this case odd 
terms have not to be computed. This difficulty can be alleviated, especially when 
high-speed computers are used, either by programming the differentiation ope- 
ration, or by using the Runge-Kutta method. Differentiation has already been 
programmed by KAHRIMANIAN [4], in the form of obtaining the m-th derivative 
of a function when the derivative of the elementary functions on which it depends 
are known. : 

KAHRIMANIAN makes use of Leibnitz formula of differentiation of a product 
and of Faadi Bruno formula of differentiation of a function of a function. Al- 
though this program is in itself very interesting, it seems to us that it is better _ 
to use the straightforward process of determining the sucessive derivatives one ~ 
from the other, even if this will necessitate the combination of similar terms at 
the end of the computations. 

The method of RuncE Kurta, consists in computing in succession 


ky=h ft, %9) 
=h f(t) + ah, % + BR), 
hg = h fltg + a hy %) + By ky + By Re) , 


and to determine £, — % as a linear combination of the h’s. It turns out that four 


o 
| 


Vol. VIII, 1957 Kurze Mitteilungen — Brief Reports — Communications bréves 155 


computations give 11 equations in 13 unknowns which have a double infinity 
of solutions. In a report [2] we discuss this méthod to start the computation in 
detail; we also use the Richardson metamethod to obtain greater precision as well 
as a check; finally we program the complete method, in such a way that it can 
be given to a coder with little mathematical experience whatever be the digital 


computer used. 
5. On a Second Order Differential Equation 


As we have mentioned already in 3 (i), a treatment similar to the above one 
can be made for many other difference methods. We will, just to provide an 
exercise, state the following theorem: 


Theorem III. /f f(t, «) has continuous n-th partial derivatives in a domain D, 
the solution x of the differential equation 


= f(t, x) (5.1) 


with initial conditions (ty, %))€ D and (ty, %y) is approximated by the solution of 
the difference equation. 


O*S Sih {bee ct tyke) = 495 (5. 2) 
z ) Note Bonk tat 
f= %o t+ hxyt+ = %+ Peed Nom ret Wed a oes Brey a eet eo | 
S 1 (5, 3) 
Jimeno. | 
in such a way that 
E(t) = #(t) + A my(t) + ++ + A? y(t) + ANF? CE, h) , (5.4) 


where 7,(¢) ave continuous functions and €(t;, h) tends to zero with h. In other words 
the difference method is of class RC(n + 1), (see [1]). 
More specifically 7,(¢) is the solution of 


: 0 1 ; 
Pes — UF ae #'°, (to) = tilts) = 9, (5. 5) 
- of 1 0? Leaky) 1 F ne rng oy 2 
Ne ag Me 2 x ey ni 12 ees 360 #°% , Molto) = Neto) = 9. (5. 6) 
Joun Topp does give [6] an example of the use of (5.2) to solve ¥= —4, 


x(0) = 0, #(0) =1. His five decimal places computations [Table I, (5)] do start 
with the exact value of ¥ at ¢= 0-1, instead of the value (5.3) which should be 
taken. The error is here negligeable because +*?(0) = 0 and 


1 oe La 
ete ete ae = 14% 10-7! 
(0-1)8 (s+ ia) x°(0) 
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Résumé 


Dans cet article nous attirons l’attention sur un résultat di a GAUNT et que 
voici: 
Si pour résoudre |’équation différentielle 


& = f(t, *), 4(0) = % 
on utilise la méthode des différences 
weg; =hflt, €)¢, = %o & 


la valeur de é, doit étre calculée d’une fagon déterminée, distincte de *#(h). 
L’erreur varie alors uniformément au sens décrit ci-dessus et on peut appliquer 
la métaméthode de RICHARDSON. 

Nous donnons un énoncé précis et une démonstration élégante que nous 
étendons ensuite au cas d’une équation différentielle de second ordre. Nous 
indiquons d’autres généralisations possibles et signalons une méthode pratique 
de calcul pour &,. 


(Received: August 20, 1956.) 


Jet Issuing in All Directions From a Thin Slit Round a Circular 
Cylinder of Small Radius Containing Liquid Under Pressure 


By LaksHMI SANYAL, Calcutta?) 


Two infinite circular cylinders of equal cross-sections are placed very near 
each other with their axes along the axis of zg and contain liquid under pressure. 
The liquid issues through the opening in all directions at right angles to the axis, 
and mixes with the surrounding liquid of the same kind at rest. We take the 
origin on the axis at the level of the slit, and use cylindrical coordinates 7, @, z. 
The motion in the jet is symmetrical about the z-axis. 

Let u, 0, w be the components of velocity in the directions of 7, 0, z.uandw . 
are independent of 0. Assuming the motion to be steady, Stokes-Navier equations 
in cylindrical coordinates reduce to two, viz., 


Ou Ou il Op 0?u di Ou 0?u u 
ae + 5 =X ae th ge | GK Gee eee oe 
Ow Ow 1 Op Wo, a dw 0?w 
Oy Seats Oe Q Opa ( Ove’ > “By By8 ) (2) 


1) University College of Science. 
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The equation of continuity is 


0 0 

= (yu) + ae (yw)=0. (3) 

Assuming u, 0u/Or, 0?u/dv? to be of o(1) in the jet, we see from the equation 

of continuity that 0(y w)/0z = 0(1); therefore w = 0(6), where 0 is the thickness 

of the jet. Since w= o0(1), du/dz = 0(1/6) within the jet. The left-hand side of 
equation (1) is therefore of o(1). On the right-hand side the terms 


O2u 1 Ou u 


Ore? ox Oe? 4h 
are of 0(1), while 0?%/dz* is of 0(1/62). Therefore the coefficient of » on the right- 


hand side of equation (1) is of 0(1/d*). If the motion in the jet is to be influenced 
by viscosity, »y must be of 0(d?). Hence the first equation of motion reduces to 


Ou. Ou 1 Op 07u 


ieee ee Ss ae Saat 


If we consider now the second equation of motion, its left-hand side is seen 
to be of 0(6). On the right-hand side 


0?w xt Ow 
“Or 7 oy Or 


are of 0(0), while 0?w/0z? is of 0(1/6), so that 


- 


es 1 oe. 02w 
*\ Or? Y x Or ar) 


is of o(0). 
Hence we get from the second equation of motion 
op 
rs eee 


If we neglect terms of 0(6), the pressure in the jet can be considered to be a 
function of y alone and for a given value of 7, it is given by the pressure of the 
liquid outside the jet. Since we consider the liquid outside the jet to be at rest, 
we have 0p/dv = 0. Therefore the equation of motion in the jet becomes 

0 Ou 07u 


uU 
— : 4 
ar ae eae ” Oz? re) 


If y be the stream function of the flow in the jet 


1 Ow i ow 
le ilo : 5 
Ss Y Oz” - Y x Gy ) 


To obtain a similarity solution of equation (4), we assume 


Pe a oe (6) 


yt 
where a, b, p, g are constants. Then 


u=abr?-I-l1f, w=—ar?-*(pf— qf). 
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Substituting in equation (4) and simplifying, we get 
BYP HIRD ig tea ep fale ie (7) 
In order that this equation may reduce to an ordinary differential equation 


in f and € only, we must have 
: 2 Ded Sees (8) 


The momentum which crosses a section of the jet perpendicular to v = const, 
between two axial planes inclined at an angle d@ is 


M a0 = 297d0 | ude 2 gab a0 re>-a-+ | fae. 
0 


Since the motion is steady and the pressure is constant, /d@ must be indepen- 
dent of v. Therefore 


2p G1 aie (9) 
From (8) and (9), we have 
Substituting in (7) and taking a = 2 yb, we see that f satisfies the equation 
2+ ff) + f= 0. (11) 


From symmetry, y = constant and du/0dz = 0 when z= 0,1. e., 
{Oy = OF OP(O) Os 
Integrating (11) with these conditions, we get 
(rat a 


This equation is identical with the corresponding equation obtained by BIcKLEY 
in his solution of the problem of the two-dimensional jet. Integrating again, we 


get 
j= 7 const = 62 


f=ctanh (¢€+ B). 


Since /(0) = 0 we have f = 0. The condition u > 0 when z > ov, i.e., f’(£) > 0 
when € + oo, is satisfied. We have therefore 


which gives 


= acy tanh(~2*). 


Replacing ac, bc by a,b respectively we do not alter the value of the ratio 
b/a = 1/2 y. This means that we can take c = 1. We then have 


be 


/(¢) = tanh¢t ~ p=artanhe, ¢= =, 


RE eet ho 


Therefore if M is prescribed 
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? 


3” M \1/3 3M \1/3 z 
Dae? (—— #tanhé; CS ESA as - 
¥ ( 20 ) g ; \16 or? Y 
9 M2 \1/3 1 ‘ 
uH = (55 >) — sech2 » . mW = 3 yt M ibe 1 
32 97 ¥ 20 


The total flux of mass across the surface of the cylinder v = const, is 


my (¢ sech?¢ — tanh@) . 


oO 
Q=4nor/udz=2n (12 ety Myr, 
0 


In conclusion, I wish to express my thanks to Dr. S. Guosu for helpful suggestions. 


Zusammenfassung 


Das Problem eines zweidimensionalen Strahles, der aus einem diinnen Schlitz 
herausstrémt, ist gelést worden durch ScHLICHTING?) in Form von unendlichen 
Reihen, und durch BicKLEey’) in geschlossener Form. Das Problem eines drei- 
dimensionalen, achsensymmetrischen Strahles ist gelést worden durch ScuHLicH- 
TING?). In vorliegender Publikation wird das Problem eines dreidimensionalen 
Strahles aus einem kreisrunden Schlitz in der Wand eines Kreisrohres behandelt. 
Zu diesem Zweck wurde die Annahme getroffen, dass die Grenzschichtapproxi- 
mationen im Strahl giiltig sind. Bei Einfiihrung der Stokesschen Stromfunktion y 
ist eine Ahnlichkeitslésung der Differentialgleichung fiir y gefunden und das 
Problem in geschlossener Form gelést worden. Der totale Mengenfluss durch 
einen zylindrischen Schnitt des Strahles wurde berechnet. 


(Received: August 15, 1956.) 


Statistische Erwagungen 
bei Versuchen zur Beeinflussung der Niederschlage 


Von Pavut ScuHmID, Zirich?) 


Seit einigen Jahren wird besonders in den niederschlagsarmen Gegenden der 
Vereinigten Staaten versucht, durch Impfung der Atmosphare mit geeigneten 
Stoffen die Wolken zum Ausregnen zu bringen. An einzelnen Haufenwolken, die 
vom Flugzeug aus geimpft wurden, konnte gezeigt werden, dass unter giinstigen 
Bedingungen eine Wirkung erzielt werden kann. Man weiss aber nicht, ob die 
grossangelegten Impfaktionen, bei denen durch am Boden aufgestellte Genera- 
toren verhaltnismassig grosse Mengen von Silberjodid zerstaubt werden, eine 
ins Gewicht fallende Wirkung haben. Die einzige Méglichkeit zur Beantwortung 
dieser Frage besteht in der statistischen Auswertung der Resultate von solchen 
Impfaktionen. Zu diesem Zweck sind verschiedene Auswertemethoden angewandt 


worden. 


2) H. Scuricutine, Z. angew. Math. Mech, 13, 260 (1933). 
3) W. G. Bicxiry, Phil. Mag. 23, 727 (1937). 
1) Eidgendssische Anstalt fiir das forstliche Versuchswesen an der ETH. 
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Die einfachste Méglichkeit besteht wohl darin, dass man die Niederschlags- 
menge der Jahre, in denen ein Teil der Frontdurchgange geimpft wurde, mit der 
mittleren Niederschlagsmenge friiherer Jahre vergleicht. Da die natiirlichen Nie- 
derschlagsmengen von Jahr zu Jahr ziemlich stark schwanken, wird man mit 
diesem einfachen Verfahren aber nie zu einem Entscheid gelangen koénnen. 


Die Impfunternehmer verwenden vor allem eine sehr anschauliche Auswerte- 
methode. Uber das Testgebiet, in welchem der Niederschlag vermehrt werden 
sollte, und iiber seine weitere Umgebung sind RegenmeBstationen verteilt, an 
denen wahrend vieler Jahre periodisch die Niederschlagsmenge gemessen wird. 
Der Niederschlag eines geimpften Frontdurchganges wird ebenfalls an allen 
Stationen gemessen und daraus fiir jeden Ort der Quotient 


Niederschlag des geimpften Frontdurchganges 
Mittlere Niederschlagsmenge ohne Impfung _ 


R 


berechnet. Man nimmt an, dass die Verteilung der Niederschlage tiber das Gebiet 
ohne Impfung bei allen Frontdurchgangen ungefahr dieselbe sei. R sollte daher, 
wenn die Impfung nicht gewirkt hat, tiberall etwa gleich gross sein. Wenn die 
Impfung gewirkt hat, wird R im Testgebiet grdsser sein als in einer weiteren 
Umgebung, die nicht mehr von der Impfung beeinflusst ist; man wird dort also 
ein Maximum erhalten. Beobachtet man anderseits bei einer Impfung ein aus- 
gepragtes Maximum im Testgebiet, so wird man vermuten, dass die Impfung 
gewirkt hat. 


Abbildung 1 
Santa Barbara County, Impfergebnis fiir einen Einzelfrontdurchgang vom 10. bis 13. Januar 1952, 
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weave a 


Trusreremsiurn = Reece? 
Precipitation Gages: * 
Stationary generators & 


Long record stations: => 


Abbildung 2 
Santa Barbara County, Impfergebnis fiir die Periode vom Dezember 1952 bis April 1953. 


Die Annahme, dass die Niederschlagsverteilung ohne Impfung bei jedem 
Frontdurchgang etwa die gleiche w4re, ist nicht begriindet. Daher gibt eine solche 
Auswertung niemals ein eindeutiges Resultat. Auch ein sehr ausgepragtes 
Maximum, wie es im Beispiel der Abbildung 1?) auftritt, konnte ohne Impfung 
auftreten. i 

Oft wird die Niederschlagsmenge einer ganzen Periode, wahrend der geimpft 
wurde, als Zahler von FR beniitzt (Abbildung 2). Wenn die Impfung nicht wirkt, 
kann jetzt mit etwas grd6sserem Recht erwartet werden, dass F tiberall etwa gleich 
gross sein werde, da sich verschiedene Frontdurchgange tiberlagern. Grundsatz- 
lich kann man aber auch so nicht zu einem Entscheid gelangen. 


Bedeutend mehr diirfen wir von einer Auswertemethode erwarten, bei der die 
Regenmenge im Testgebiet mit derjenigen in einem Kontrollgebiet verglichen 
wird. Das Kontrollgebiet muss so gewahlt werden, dass es von der Impfung voll- 
kommen unbeeinflusst bleibt. Die Werte des Niederschlags im Testgebiet und im 
Kontrollgebiet werden, fiir viele impffreie Perioden, in einem Diagramm einge- 
tragen. Wenn die beiden Gebiete nahe beieinander und topographisch ahnlich 
gelegen sind, werden die Niederschlagsmengen gut korreliert sein, die Punkte des 
Diagramms werden sich durch eine Gerade ausgleichen lassen. Die beste Gerade 
wird nach der Methode der kleinsten Quadrate bestimmt. Aus der Streuung der 


2) Die Abbildungen 1-3 sind Berichten der North American Weather Consultants entnommen ; 
siehe auch Bull. S.E.V. 44, 189 (1953). 
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Punkte um diese Gerade lasst sich die Wahrscheinlichkeit 
W (y< ¥o| *); 


berechnen, dass bei gegebenem ¥ die Grésse y grosser als ein bestimmter Wert Vo ist. 
Wenn man anderseits die Griésse dieser Wahrscheinlichkeit vorgibt, zum Beispiel 
W (y<¥01*%) =1% » 


so wird durch diese Gleichung eine Funktion y,)(¥) bestimmt. Es wird im Durch- 
schnitt vieler Perioden nur etwa 1% aller Punkte oberhalb dieser Linie zu liegen 
kommen. 


E 
15 20 INCHES si 10 15 20 


Abbildung 3 
Impfgebiet Santa Barbara County. Regressionsberechnung mit zwei Vergleichsgebieten, 
Ergebnis fiir den Winter 1951/52. 
Die gestrichelten Linien sollten die Entscheidungsfunktionen (W = 1%) darstellen. In Wirklichkeit 
miissen diese Geraden aber durch nach oben gedffnete Parabelbdgen ersetzt werden. 


Wenn jetzt wahrend einer oder mehrerer Perioden geimpft wird und die 
Impfung eine ansehnliche Wirkung hat, wird der dieser Periode entsprechende 
Punkt voraussichtlich oberhalb dieser Linie zu liegen kommen. Wenn die Imp- 
fung nicht wirkt, bleibt die Wahrscheinlichkeit dafiir klein. Man hat daher grund- 
satzlich eine Moglichkeit, objektiv zu entscheiden, ob die Impfung gewirkt habe 
oder nicht. 

Diese Auswertemethode ist sehr haufig und ziemlich unkritisch angewandt 
worden. Die wahrscheinlichkeitstheoretischen Methoden, auf Grund derer die 
Entscheidungsfunktion y,(¥) bestimmt wird, beruhen auf einigen Vorausset-~ 
zungen: 

1. Der Zusammenhang zwischen x und y soll linear sein; 

2. die Verteilung von y bei festem ¥% soll normal sein; 

3. die Streuung von y soll fiir alle % gleich gross sein. 
In vielen Fallen sind diese Voraussetzungen anndhernd erfiillt. y wird um so 
besser normal verteilt sein, je grésser die Anzahl der Regenmesser ist und je 
langer die Perioden sind. Wenn die Streuung von y proportional zu ¥ anwdchst, 
werden die transformierten Variablen +* = |/x, y* = |/y die dritte Bedingung 
erfiillen. 
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Weit mehr als diese formellen Einwande failt ein prinzipieller Einwand gegen 
diese Auswertemethode ins Gewicht. Das Verhaltnis der Regenmengen in Test- 
und Kontrollgebiet wird bei jedem Frontdurchgang ein etwas anderes séin. Es 
gibt Typen von Gewittern, die dem Testgebiet verhaltnismassig mehr Niederschlag 
bringen als dem Kontrollgebiet,-auch wenn nicht geimpft wird, und solche, bei 


denen es umgekehrt ist. Wir kénnen prinzipiell drei Typen von Frontgewittern 
unterscheiden: 


Typ 4: Die Regenmenge ist im Testgebiet verhaltnismassig héher als im Kon- 
trollgebiet; 


Typ B: Die Regenmengen verhalten sich etwa so wie die Gesamtregenmengen. 


Typ ©: Die Regenmenge ist im Kontrollgebiet verhaltnismassig héher als im 
Testgebiet. 


Da das Wetter von Jahr zu Jahr und von Jahrzehnt zu Jahrzehnt seine 
Struktur etwas andert, kann es sein, dass in der Impfperiode mehr Gewitter vom 
Typ A vorkommen als in den friiheren Perioden, wodurch eine Impfwirkung 
vorgetauscht wird; es kann aber auch das Gegenteil eintreten. 

Noch schwerwiegender wird dieser Einwand, wenn nur ein Teil der Front- 
durchgange geimpft wird und nur bei diesen Gewittern die Regenmenge in den 
beiden Gebieten gemessen wird. Die Gewitter, die fiir die Impfung giinstig sind, 
k6nnen ja bevorzugt einem der drei Gewittertypen angehéren. Der Impfunter- 
nehmer kann zudem, bewusst oder unbewusst, verhaltnismassig mehr Gewitter 
vom Typ 4 impfen, wodurch natiirlich das Versuchsresultat sehr stark verfalscht 
werden kann. 

Man kénnte versuchen, nach meteorologischen Gesichtspunkten jeden Sturm 
in eine der drei Kategorien einzureihen und fiir jede Kategorie eine besondere 
Regressionsgerade aufzustellen. Das Resultat kann so etwas verbessert werden, 
die Zuordnung ist jedoch héchst unsicher und kann ebenfalls subjektiv beein- 
flusst sein. 


Wir kénnen einzig und allein dadurch zu einer objektiven Beantwortung 
unserer Versuchsfrage kommen, indem der Entscheid, ob ein Frontdurchgang 
geimpft werden soll oder nicht, dem Zufall iiberlassen wird. 

Nach meteorologischen Gesichtspunkten wird vorerst wie bisher entschieden, 
ob man das kommende Gewitter impfen sollte oder nicht. Dann entscheidet ein 
Zufallsexperiment, zum Beispiel ein Wiirfel, ob das Gewitter geimpft wird oder 
ob es als Vergleichsgewitter dient. Es werden also nur etwa halb so viele Front- 
durchgange geimpft wie bisher. Fiir jedes Gewitter werden die Regenmengen in 
den beiden Gebieten gemessen, wodurch wir zwei vollkommen gleichberechtigte 
Serien von Wertepaaren (7, y) und zwei Regressionsgeraden erhalten, von den 
geimpften und von den ungeimpften Gewittern. 

Man kann jetzt angeben, wie gross die Wahrscheinlichkeit dafiir ist, dass 
zwischen den beiden Regressionsgeraden ein bestimmter Unterschied besteht, 
unter der Hypothese, dass die Impfung keine Wirkung habe. Wenn diese Wahr- 
scheinlichkeit fiir den wirklich eingetretenen Unterschied zwischen den beiden 
Geraden geniigend klein ist, zum Beispiel kleiner als 1%, wird man daraus 
schliessen, dass die Hypothese falsch ist, dass die Impfung also gewirkt hat. 
Die Grésse der Wirkung kann einigermassen abgeschatzt werden. F 

Durch diese Versuchsanlage erreichen wir, dass Daten aus friiheren Perioden, 
in denen das Wetter ja eine etwas andere Struktur haben konnte, nicht benotigt 
werden; der Versuch ist in sich abgeschlossen. Zudem sind jetzt subjektive Ein- 
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fliisse auf das Versuchsresultat sozusagen ausgeschlossen. Bei geniigend langer 
Versuchsdauer wird daher eine dusserst zuverlassige Aussage moglich sein. 
Wenn kein passendes Vergleichsgebiet zur Verfiigung steht, kann der Versuch 
trotzdem auf diese Art durchgefiihrt werden. In diesem Fall miissen die Mittel- 
werte der beiden Vergleichsserien miteinander verglichen werden, statt der 
Regressionsgeraden. Der Versuch wird allerdings etwas langer dauern. 


Bis heute ist:einzig durch die Universitat New York ein kleiner Versuch auf 
ahnlicher Basis durchgefiihrt und abgeschlossen worden. Die Wirkung der Imp- 
fung wurde dabei aber nicht gesichert. In Israel lauft seit einigen Jahren ein 
Versuch, bei dem durch ein Zufallsexperiment bestimmt wird, welche Gewitter 
geimpft werden und welche als Vergleichsgewitter dienen sollen. Das Zufallsexpe- 
riment wird aber zum voraus durchgefiihrt, also vor dem Entscheid des Meteoro- 
logen. Da der Meteorologe weiss, ob das nachste von ihm bestimmte Gewitter 
geimpft wird oder nicht, kann er das Versuchsresultat subjektiv beeinflussen, 
indem er zum Beispiel wartet, bis ein Gewitter vom Typ 4 eintritt. 

Wir haben daher bis heute noch keine Gewissheit dariiber, ob die getatigten 
Impfaktionen eine messbare Wirkung haben. Auf Grund der bisher durchge- 
fiihrten Auswertungen werden mittlere Wirkungen von 6-15% der mittleren 
Regenmenge vermutet. Diese Differenz kann aber sehr wohl durch subjektive 
Einfliisse vorgetauscht werden. 


Die gleichen Impfmethoden werden angewendet zur Bekampfung des Hagels. 
Grundsatzlich kénnen in diesem Fall genau die gleichen Uberlegungen gemacht 
werden. Auch hier wird man nur zu einem eindeutigen Entscheid gelangen kénnen, 
wenn ein Zufallsexperiment bestimmt, ob ein Gewitter geimpft werden soll oder 
nicht. Ein Kontrollgebiet, das unbeeinflusst von der Impfung ist, dessen Hagel- 
falle aber gut korreliert sind mit denjenigen im Testgebiet, ist meist nicht vorhan- 
den, so dass man sich darauf beschranken muss, die totale Hagelhaufigkeit in den 
beiden Gewitterserien zu vergleichen. Eine weitere Erschwerung entsteht dadurch, 
dass die Hagelhaufigkeit oft sehr klein ist. Immerhin darf dafiir vielleicht eine 
etwas hohere prozentuale Wirkung der Impfung erwartet werden, was die Ver- 
suchsdauer betrachtlich verkiirzen kénnte. Zuverlassige Resultate liegen bis 
heute aber auch tiber diese Impfwirkung noch nicht vor. 


Zusatz bei der Drucklegung: Soeben wird bekannt, dass im Santa Barbara 
County, auf das sich auch die obigen Figuren beziehen, seit 1. Januar 1957 ein 
Regenvermehrungsversuch mit Zufallszuordnung im Sinne dieser Ausfiihrungen 
lauft. Die Impfung wird wie bis anhin durch die North American Weather Con- 
sultants besorgt, und die Niederschlagsmessung untersteht dem Amt der California 
State Water Resources. Die Versuchsplanung, die statistische Auswertung sowie 
das Zufallsexperiment ist vom Statistical Laboratory der University of California, 
Berkeley, unter Leitung von J. NEyman, tibernommen worden. Auch in dem von 
der Eidgendssischen Kommission zum Studium der Hagelbildung und Hagel- 
‘abwehr im Siidtessin durchgefiihrten Impfversuch zur Behinderung der Hagel- 
bildung wird diesen Sommer das gleiche statistische Verfahren angewandt werden. 
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Summary 


This report deals with the experiments toincrease precipitation and to 
prevent hail by way of seeding the atmosphere with chemicals suited to the 
purpose. The statistical evaluation methods hitherto adopted are discussed, and 
it is shown why it is impossible to give a reliable answer to the question of the 
seeding results unless randomization is introduced into the scheme of the experi- 
ments. 


(Eingegangen: 12. Januar 1957.) 


The Transfer of Momentum Across Surfaces Drawn Perpendicular 
to the Direction of Flow in Gases Under Laminar Conditions 


By IBRAHIM IBRAHIM SHERIF, Alexandria, Egypt?) 


It is not generally remarked in treatises on kinetic theory of gases that 
tangential stress forces act across surfaces that are drawn perpendicular to the 
direction of flow in gases under laminar conditions. The purpose of the present 
work is to calculate these forces and to prove that they must be equal to the 
stress forces that act across surfaces parallel to the direction of flow; otherwise 
a.cube of fluid with faces perpendicular to the two directions would obviously 
be set into rotation. 

The transverse momentum across surfaces perpendicular to the direction of. 
flow arises in the following way: 

Suppose the flow is horizontal and increases in the upward direction. Consider 
a face of an elementary cube across which the mass flow is perpendicularly 
outward from the element. The number of molecules entering the element per 
second across this face due to thermal agitation is slightly decreased by the mass 
velocity. This decrease in the number is greater in the case of molecules that enter 
with a downward component of velocity than in the case of molecules that enter 
with an upward component, because the former molecules come from a region of 
the gas in which the mass velocity is greater. Thus molecules entering with an 
upward component are slightly more numerous, and the result is a net inflow 
of upward momentum into the element. This effect is easily calculated as follows: 

Let us first calculate the shearing stress across surfaces parallel to the direction 
of flow: 

Let N be the number of molecules per cubic centimeter, therefore number of 
molecules per unit volume with thermal velocity C in the direction d0 dp is 


sin 0 
4a 


N do dq; 


component of C for these molecules in y-direction is 
C sin 6 cos 9; 
the stream velocity contributes nothing in the y-direction, hence number of 


1) Faculty of Engineering, University of Alexandria. 
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s ae 
molecules crossing unit area per second in the y-direction is 


we 


N eae d0 dy C sin6 cos; 


momentum in #-direction for one molecule is 


du. 
| C cos@ + (u- Ls sin 9 cos ¢| \, 


OS 


Figure illustrating momentum. 


where wu is the stream velocity and L is the mean free path. Hence total momen- 
tum transfer is: 


Nm | | Csind cos (Ccosd + w— L ie sind cos) ee d0 dp 
* (1) 
=-NmCZL ie sin? d9 = — > eCL 
0 


Second, the shearing stress across surfaces perpendicular to the direction of flow 
can be calculated in a similar manner as follows: 
Component of C for the molecules in the v-direction is: 


Ccosé, 


these molecules near the element have also a velocity in the x-direction of average 
value 


du. 
US Bey, sin cos , 


- 


an - Ne . a “ : * uy i 
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_ hence number crossing unit area per second in the x-direction is 
oa 


te ” 


(c cosé+u—L a sin 6 cos 7) N a dadp, 


momentum in y-direction is: “ 
; mcsin@ cosy , 


the stream component contributing nothing, hence momentum transfer is 


2% 7 
N ey ewer mit Saeatny ete sin @ 
m | / C sin @ cos p (Cc cos@+u—L ae sin 6 cos p rie d0 dp 
0 0 (2) 
1 NS GU 


This means that the momentum across surfaces parallel to the direction of flow 
is equal to the momentum across surfaces in the perpendicular direction. In 
other words the principle of conjugate shears which is generally supposed to 
apply to elastic solids must as well hold in gases. 


Zusammenfassung 


Bei laminarer Str6mung eines Gases ist die Ubertragung der Bewegungsgrisse 
quer zu Flachen senkrecht zur Str6mungsrichtung berechnet und gleich der 
transversalen Bewegungsgrésse quer zu Flachen, die parallel zur Stro6mungs- 
richtung laufen, gefunden worden. Daraus kann der Schluss gezogen werden, 
dass ein Fliissigkeitswiirfel mit Seiten, die parallel zu den angegebenen Richtungen 
sind, statisch im Dreh-Gleichgewicht ist. 


(Received: September 5, 1956.) 
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Properties of Metallic Surfaces. A Symposium Organized by the Institute 
of Metals and Held at the Royal Institution, London, on 19 November 1952 ([Insti- 
tute of Metals Monograph and Report Series No 13] The Institute of Metals, 
London 1953). 368 S., 150 Abb.; 35s. 

In diesem Symposium wurden von hervorragenden englischen Wissenschaf- 
tern (BowpEN, ToLANSKy, TaBor, Evans und anderen) neue Ergebnisse von 
Grundlagenforschungen im Gebiet der Oberflachentextur und -struktur von 
Metallen bekanntgegeben, die vielfach einen wertvollen Beitrag zur Klarung 
physikalischer, chemischer und technologischer Eigenschaften von Metallen bil- 
den. Diese Grundlagenforschung ist nicht nur fiir Physiker und Chemiker von 
héchstem Reiz, sondern erweckt auch das lebhafteste Interesse des Technikers, 
der sich bewusst ist, wie ausserordentlich problematisch solche Rechnungsbei- 
werte sind, wie Reibungskoeffizienten und andere mehr. 

Die Veroffentlichung der Vortrage und Diskussionsbeitrage von gut zwei 
Dutzend englischen Forschern verdient deshalb héchste Beachtung aller, die 
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aus physikalisclf® theoretischem oder technisch-praktischem Interesse sich mit 
den Vorgangen ‘befassen, die von den dussersten Grenzschichten ausgehen oder 
durch sie stark beeinflusst werden, seien es Probleme der Reibung, der Korrosion, 


der Wechselfestigkeit, der Diffusion, des Skineffektes, der Beilby-Schicht oder — 


andere mehr. 
Die Hunderte von Literaturangaben geben ein eindriickliches Bild von dem 


Umfang, den diese «Wissenschaft der Grenzschichten fester K6rper» heute ange- 
nommen hat, deren Auswirkung im Gebiet der praktischen Technik nur noch 
eine Frage der Zeit ist. E. Bickel 


Mass und Integral und ihre Algebraisierung. Von C. CARATHEODORY 
(Birkhauser Verlag, Basel und Stuttgart 1956). 337 S., 12 Fig.; Fr. 38.50. 
Die den modernen Integraltheorien zugrunde gelegten Punktmengen haben 


als wesentliche Eigenschaft, dass sie als Operationen die Bildung von Vereinigung, 


Durchschnitt und Differenz zulassen. Dieselben Eigenschaften weisen nun aber 
auch die Elemente (Somen) einer Booleschen Algebra auf. Geht man daher von 


einer solchen Algebra aus, dann erhalt man eine die bisherigen Integraltheorien 


umfassende Theorie von grosser Allgemeinheit und Natiirlichkeit, und man 


gelangt zu einer Algebraisierung des Integralbegriffs. Diese Zusammenhange | 


wurden von C. CARATHEODORY vor ungefahr zwanzig Jahren bemerkt. — Das 
Manuskript zum vorliegenden Buch, das diese Ideen tibersichtlich entwickelt, 
konnte CARATHEODORY noch vor seinem Tode fertigstellen. Die Bereinigung fiir 
den Druck lag in den bewahrten Handen von P. FINSLER, A. ROSENTHAL und 
R. STEUERWALD, denen es zu verdanken ist, dass das Buch tiberhaupt erscheinen 
konnte. 

Am Anfang wird eine Boolesche Algebra eingeftihrt, unter Zugrundelegung 
des Axiomensystems von M. H. Stone, und dann werden die Begriffe Ortsfunk- 


tion, Massfunktion und Integral definiert und ihre Eigenschaften entwickelt. 


Schliesslich folgen Anwendungen der Integraltheorie insbesondere auf Grenz- 
prozesse, Mass- und Inhaltsfunktionen. — Die Darstellung ist ausgezeichnet und 
tibersichtlich, und das Buch gibt die erste treffliche lehrbuchmassige Einfitihrung 


in das Gebiet. Ebenso sind auch Druck und Ausstattung des Buches in gewohnter - 


Weise vorziiglich. E. Roth-Desmeules 


Vorlesungen uber Himmelsmechanik. Von C. L. SircEL ({Grundlehren © 


der mathematischen Wissenschaften, Bd. 85] Springer-Verlag, Berlin, Géttingen, 
Heidelberg 1956). 212 S.; DM 29.80/DM 33.-. 

Das vorliegende Werk von C. L. SIEGEL behandelt einige Ideen und Resultate 
der klassischen Himmelsmechanik, die aber im Grunde nicht nur dort ihre 
Bedeutung haben, sondern sich auch auf andere Falle iibertragen lassen und 


vielfach zu wenig bekannt sind. Zunachst geht es um das Verhalten der Lésungen 
von gewohnlichen Differentialgleichungen im Grossen; hier haben H. PoINcARE — 


und K. F. SunpMAN sowie G. D. BIRKHOFF wesentliche Beitrage geleistet. Ein 
Kapitel befasst sich mit der Bestimmung von geschlossenen (periodischen) Lésun- 
gen, ein anderes mit dem Stabilitatsproblem (Satz von LyapuNov) und der damit 
zusammenhangenden, so schwierigen und weitgehend noch ungelésten Frage der 
Konvergenz der auftretenden Entwicklungen. 

Die Darstellung des Verfassers ist sehr klar und iibersichtlich, und das aus 
Vorlesungen hervorgegangene ausgezeichnete Buch vermag jedenfalls auch 
manche Anregung zu analogen Untersuchungen auf andern Gebieten als der 
Himmelsmechanik zu geben. E. Roth-Desmeules 
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Neuerscheinung 


RUTISHAUSER 


Der Quotienten-Differenzen-Algorithmus 
Von Prof. Dr. H. Rutishauser, ETH, Ziirich. 74 Seiten. Broschiert Fr.8.50 (DM 8.50). 
Mitteilungen a. d. Institut fiir angew. Math. ETH. — Nr.7. 


Der Quotienten-Differenzen-Algorithmus ist eine neu entwickelte Methode zur Lésung 
verschiedenartiger Aufgaben der angewandten Mathematik, von denen nur die 
folgenden genannt seien: Summation schlecht konvergenter Reihen, Lésung alge- _ 
braischer Gleichungen, Bestimmung von Eigenwerten vielreihiger Matrizen, insbe- 
sondere mit Rechenautomaten. Er hat sich fiir die Bestimmung kritischer Frequenzen 
und Drehzahlen bewahrt. Der Verfasser gibt eine zusammenfassende Darstellung. 


Voranzeige 


WITTENBERG 


Vom Denken in Begriffen 


Mathematik als Experiment des reinen Denkens. Von Alexander Israel Wittenberg, 
Professor an der Laval-Universitat in Quebec. Etwa 370 Seiten. Preis gebunden 
ca. Fr. 25.— (ca. DM 25.-). Sammlung « Wissenschaft und Kultur» — Band 12. 


Nirgends wird die Macht und die innere Koharenz begrifflichen Denkens so scho- 
nungslos auf die Probe gestellt wie in der Mathematik. Im mathematischen Grund- 
lagenproblem wird, wie dieses Werk in eingehender Diskussion zeigt, eine echte 
Problematik menschlichen Denkens offenbar. Die Untersuchung derselben fiihrt zu 
wesentlichen Einsichten in die Natur und die Grenzen menschlicher Erkenntnis, 
und damit iiber den Menschen schlechthin, 
Das Buch will neues Licht auf die mathematische Grundlagenforschung werfen. 
Zugleich wendet es sich an alle, insbesondere Wissenschaftler, Philosophen und 
Theologen, die sich tiber die erkenntnismassigen Gegebenheiten des menschlichen 
Daseins Rechenschaft ablegen wollen. 
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BIRKHAUSER VERLAG: BASEL UND STUTTGART 


v. MANGOLDT-KNOPP 


Einfihrung 
in die hohere Mathematik 


10. vollstindig neubearbeitete Auflage 


BAND I: Zahlen, Funktionen, Grenzwerte, Analytische Geo- ~ 
metrie, Algebra, Mengenlehre. XVI + 564 Seiten mit 116 
Figuren, Gr. 8°. Ganzleinen DM 23.-. 1956. 


BAND II: Differentialrechnung, Unendliche Reihen, Ele- 
mente der Differentialgeometrie und der Funktionstheorie. 
XVI + 630 Seiten mit 108 Figuren. Ganzleinen DM 23.-. 
Februar 1957. f 


BAND III: Integralrechnung und ihre Anwendungen, -Funk- 
tionentheorie, Differentialgleichungen. XVI + 618 Seiten mit 
103 Figuren. Gr. 8°. Ganzleinen DM 23.—. Juli 1957. 


Niederfrequenz- und 
Mittelfrequenz-Messtechnik 


Von Dr.-Ing. A.Wirk und Dipl!-Ing. H. G. THILo. 
VIII + 234.Seiten mit 223 Abbildungen. Gr. 8°. 
Ganzleinen DM 28.-. 1956. 


Das Buch gibt ein vollstandiges Bild der Physik und Technik 
dieser Messverfahren und behandelt auch Entwicklung, Bau, Prii- 
fung, Montage und Wartung solcher Anlagen und Messgerate. 


S.HIRZEL VERLAG STUTTGART 


